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ABSTRACT 

We present results of a spectroscopic search for Lya emitters (LAEs) in the C11604 supercluster field 
using the extensive spectroscopic Keck/DEep Imaging Multi-Object Spectrograph database taken as 
part of the Observations of Redshift Evolution in Large Scale Environments (ORELSE) survey. A 
total of 12 slitmasks were observed and inspected in the C11604 field, spanning a survey volume of 
1.365 X lO'' co-moving Mpc^. We find a total of 17 high redshift (4.39 < z < 5.67) LAE candidates 
down to a limiting flux of 1.9 x 10~^® ergs s~^ cm^^ (L(LyQ!) = 4.6 x 10''^ ergs s~^ or ^^O.IL* at 
z ^ 5), 13 of which we classify as high quality. The resulting LAE number density is nearly double 
that of LAEs found in the Subaru deep field at z ^ 4.9 and nearly an order of magnitude higher 
than in other surveys of LAEs at similar redshifts, an excess that is essentially independent of LAE 
luminosity. We also report on the discovery of two possible LAE group structures at z ~ 4.4 and 
z ~ 4.8 and investigate the effects of cosmic variance of LAEs on our results. Fitting a simple 
truncated single Gaussian model to a composite spectrum of the 13 high quality LAE candidates, we 
find a best-fit stellar velocity dispersion of 136 km s~^. Additionally, we see modest evidence of a 
second peak in the composite spectrum, possibly caused by galactic outflows, offset from the main 
velocity centroid of the LAE population by ~440 km s~^ as well as evidence for a non-trivial Lya 
escape fraction. We find an average star formation rate density (SFRD) of ~ 5 x 10^'^ M© yr~^ 
Mpc^'^ with moderate evidence for negative evolution in the SFRD from z ~ 4.6 to z 5.7. By 
simulating the statistical flux-loss due to our observational setup we measure a best-fit luminosity 
function characterized by = 2.2^^t^'^^ x 10^^ ergs s~^ Mpc^^ for a=-1.6, generally consistent with 

measurements from other surveys at similar epochs. Finally, we investigate any possible effects from 
weak or strong gravitational lensing induced by the foreground supercluster, finding that our LAE 
candidates are minimally affected by lensing processes. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: clusters: general — galaxies: 
high-redshift — techniques: spectroscopic 



1. INTRODUCTION 

While Lya emitters (LAEs) have been sought for 
nearly 40 years, designing and implementing surveys ca- 
pable of detecting large unbiased populations of these 
objects have proven difficult. Due to the extreme faint- 
ness of the population and technological limitations, 
the searches pioneered by Davis et al. in the 1970s 
(Davis & Wilkinson 1974; Partridge 1974) established 
what would later be a theme for such surveys: con- 
straints on galaxy populations and cosmological param- 
eters through a dearth of detections. At that time little 
was known about the properties of high-redshift galax- 
ies, with the observational distinction between LAEs and 
a second high-redshift star-forming population, Lyman 
break galaxies (LBGs), not yet possible. This ignorance 
about the fundamental differences in the properties of 
the two types of high-redshift galaxies resulted in the 
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grouping of both galaxy populations into a single cate- 
gory: Primeval Galaxies (PGs). While early theoretical 
modeling (see Davis 1980) predicted the density of PGs 
to be > 10000 per deg^ at high redshift (z > 3), early 
searches for PGs (Koo & Kron 1980; Saulson & Boughn 
et al. 1982; Boughn et al. 1986; Pritchet & Hartwick 
1987, 1990; Elston et al. 1989; de Propris et al. 1993; 
Thompson et al. 1995; Thompson & Djorgovski 1995) 
were unable to find any such objects. It was not until 
the mid-1990's with the searches of Steidel and collabo- 
rators that large populations of PGs were detected, al- 
most exclusively of the LBG flavor (Steidel et al. 1996a, 
1996b). 

The detection of LAEs has proven significantly more 
problematic than LBGs due to the difficulty of efficiently 
identifying the Lya line in candidate galaxies. In addi- 
tion, the Lya line is only observed in ^^25% of high red- 
shift star-forming galaxies (Steidel et al. 2000; Shapley 
et al. 2003). Due to these difficulties, it is only in the 
past half-decade that techniques have been successfully 
developed and implemented to detect reasonably large 
numbers of LAEs. 

The most common technique in contemporary LAE 
searches is the use of custom-made narrowband filters 
with bandpasses of 100 A or less, designed to collect 
light in windows of low atmospheric transmission. Imag- 
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ing campaigns using such filters have been successfully 
undertaken in blank fields complemented by deep broad- 
band photometry (Hu et al. 2004, hereafter H04; Ouchi 
et al. 2003, 2008, hereafter O03, O08; Rhoads et al. 
2000; Malhotra & Rhoads 2002) or in areas of suspected 
overdensities (Kurk et al. 2004; Miley et al. 2004; Ven- 
emans et al. 2004; Zheng et al. 2006; Overzier et al. 
2008). While this technique has proven capable of de- 
tecting large numbers of LAEs, the populations detected 
may be inherently biased, due either to the small redshift 
windows probed, a bias intensified by the high level of 
observed spatial clustering of LAEs, or due to the large 
line equivalent widths (EWs) necessary to detect such 
objects. 

An alternative is dedicated spectroscopic campaigns in 
blank fields (Crampton & Lilly 1999; Martin & Sawicki 
2004, hereafter MS04; Tran et al. 2004, hereafter T04; 
Martin et al. 2008, hereafter M08), yielding samples of 
LAEs complementary to photometric searches. While 
narrowband imaging surveys provide large samples as a 
result of their ability to probe large volumes in relatively 
short periods of time, the increased sky noise due to the 
large filter bandpass (~100A) relative to a "typical" Lya 
emission width (10-20 A full-width at half-maximum, 
FWHM) makes it difficult to probe deep into the LAE 
luminosity function. As a result the line luminosities of 
galaxies detected in these surveys are usually at or above 
L*. By dispersing the night sky background so that the 
emission line has only to exceed the background over the 
natural width of the line rather than over ~ 100 A, spec- 
troscopic surveys for LAEs become much more efficient 
probes of sub-L* galaxies at high redshift. 

The difficulty with such observations is that spec- 
troscopy probes a significantly smaller area on the sky 
than narrowband techniques, with the area reduced by 
the ratio of the slit area to the telescope field of view (see 
discussion in M08) . The early dedicated searches of T04 
and MS04 suffered from this limitation, covering 17.6 
arcmin^ and 5.1 arcmin^ respectively. Along with the 
small spectral bandpasses designed to fit in atmospheric 
transmission windows, this effect severely limited the vol- 
ume probed by such surveys and as a result no LAEs were 
detected. It was not until the recent search of M08, us- 
ing similar techniques but with a significant increase in 
sensitivity and field of view, that LAEs were discovered 
exclusively through dedicated spectroscopic techniques. 
These results demonstrate the necessity of large vohime 
searches to effectively detect and analyze populations of 
LAEs. 

With the recent use of multi-object spectrographs for 
large surveys of galaxies at intermediate redshift (e.g., 
DEEP2, VVDS) it has become possible to obtain deep, 
high resolution spectra of large patches of blank sky and 
move beyond single serendipitous discoveries of LAEs 
(Franx et al. 1997; Dawson et al. 2002; Stern et al. 
2005) to statistical samples of high-redshift emission line 
galaxies (Sawicki et al. 2008; hereafter S08). With this 
in mind, we have searched the extensive (3.214 arcmin^, 
1.365x10^ Mpc^) spectroscopic database of the C11604 
superclustcr at z ^ 0.9 (Gal et al. 2008, hereafter 
G08). This structure is studied as part of the Obser- 
vations of Redshift Evolution in Large Scale Environ- 
ments (ORELSE) survey (Lubin et al. 2009), an ongo- 



ing multi-wavelength campaign mapping out the envi- 
ronmental effects on galaxy evolution in the large scale 
structures surrounding 20 known clusters at moderate 
redshift (0.6 < z < 1.3). While the angular coverage is 
moderate compared to other such surveys of LAEs, the 
C11604 data have the advantage of large spectral cover- 
age (see Section 2) and deep observations on the Keck 
10-m telescope, which allow us to probe down to un- 
precedented levels in the luminosity function (~0.1i* at 
^; ~ 5). As a result we find 17 LAE candidates in our 
moderately sized volume, almost all of which are fainter 
than the characteristic luminosity at z ^ 5. These de- 
tections allow us to place some of the first constraints 
on the properties of low luminosity galaxies at high red- 
shift, including implications for this population's role in 
the reionization of the universe. 

The remainder of the paper is organized as follows: 
Section 2 describes the spectral data and our selection 
process. Section 3 describes tests to validate our high 
redshift LAE candidates. Section 4 includes a discus- 
sion of other properties, such as photometric limits, line 
equivalent widths, velocity profiles, and star formation 
rates (SFRs) of the LAE candidates. In Section 5 we de- 
scribe the number density and luminosity function of our 
LAE candidates as well as the effects of LAE clustering 
and cosmic variance. In addition, since these data were 
taken in an area of the sky with a rare, massive structure 
in the foreground, we also discuss in Section 5 any pos- 
sible contributions from gravitational lensing. Section 6 
summarizes our results. Throughout this paper we use 
the concordance ACDM cosmology with Hq = 70 km 
s"\ r^A = 0.7, and CIm = 0.3. At z = 4.8, the median 
redshift of our sample, the age of the universe is 1.2 Gyr 
and the angular scale is 6.41 kpc arcsec"^, with 621 Myr 
elapsing between z = 6.4 and z = 4.1, the redshift range 
of LAEs to which our spectral coverage is sensitive. All 
EW measurements arc given in the rest frame and all 
magnitudes are given in the AB system (Oke & Gunn 
1983; Fukugita et al. 1996). 

2. DATA 

The first target of the ORELSE survey, and the subject 
of study in this paper, is the C11604 field, containing the 
C11604 superclustcr at z = 0.9: a massive collection of 
eight or more constituent groups and clusters spanning 
13 comoving Mpc in the transverse dimensions and 
nearly 100 comoving Mpc in the radial dimension 
(see G08 for the coordinates and velocity centroids of 
the clusters that comprise the C11604 superclustcr). The 
data on this structure include Very Large Array (B-array, 
20 cm), Spitzer IRAC (3.6/4.5/5.8/8.0 nm) and MIPS 
24 /im imaging, archival Subaru V-band imaging, deep 
Palomar r' i' z' Kg imaging, a 17 pointing Hubble Space 
Telescope ACS mosaic in F606W and F814W, and two 
deep (50 ks) Chandra pointings. 

In addition to the photometric data, an extensive spec- 
troscopic campaign has been completed in the C11604 
field to determine the rest-frame optical/UV spectral 
properties and rcdshifts of a large fraction of the con- 
stituent cluster members. Photometric data alone are 
not ideal for this purpose, as typical photometric red- 
shift errors can span the line-of-sight extent of large scale 
structures such us C11604, leading to severe uncertain- 
ties in environmental indicators such as local density. To 
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accurately quantify environmental effects, large spectro- 
scopic coverage is essential in minimizing the effects of 
projections (see G08 for a more detailed discussion). 

To this end, 12 masks covering a large portion of the 
C11604 structure were observed with the DEep Imag- 
ing Multi-Object Spectrograph (DEIMOS; Faber et al. 
2003) on the Keck II 10-m telescope between May 2003 
and June 2007. The observations were taken with 1" slits 
with the 1200 1 mm~^ grating, blazed at 7500 A, result- 
ing in a pixel scale of 0.33 A pix~^, a resolution of ^ 1.7 
A (68 km s^^), and typical wavelength coverage of 6385 
A to 9015 A. Each DEIMOS mask contained between 80 
and 130 individual slits with an average length of 9.9", 
with 95% having slit lengths between 4.92" and 14.88". 
The slits in each mask combined for a total sky coverage 
of 0.2678 arcmin^ per mask, independent of the number 
of slits. The spectroscopic targets for these slits were 
selected based on the likelihood of being a cluster mem- 
ber, determined through a series of color and magnitude 
selections (see G08). The masks were observed with dif- 
fering total integration times, which varied depending on 
weather and seeing conditions, in order to achieve simi- 
lar levels of redshift completeness of targeted galaxies. A 
differing number of 1800s exposures were stacked for each 
mask, with total integration times of 7200s to 14400s. 

The exposures for each mask were combined using 
the DEEP2 version of the spec2d package (Davis et al. 
2003)^. This package combines the individual exposures 
of the slit mosaic and performs wavelength calibration, 
cosmic ray removal and sky subtraction on slit by slit ba- 
sis, generating a processed two-dimensional spectrum for 
each slit. The spec2d pipeline also generates a processed 
one-dimensional spectrum for each slit. This extraction 
creates a one-dimensional spectrum of the target, con- 
taining the summed flux at each wavelength in an op- 
timized window. In all, 903 total high quality (Q > 3, 
see G08 for an explanation on the quality codes) spectra 
were obtained, with 329 falling within 0.84 < z < 0.96, 
the adopted redshift range of the supercluster. 

2.1. Searching for Serendipitous Detections 

During the reduction process spec2d also determines if 
any other peaks exist in the spatial profile of the slit 
that are distinct from the target. If such peaks ex- 
ist, spec2d does similar extractions at these spatial lo- 
cations creating one-dimensional spectra for these non- 
targeted serendipitous detections (hereafter serendips). 
All serendipitous spectra generated in this manner were 
systematically inspected by one of us (RG) to determine 
whether these extractions contained genuine stellar or 
galactic signatures rather than instrumental or reduction 
artifacts. 

In addition to the spec2d extraction algorithm for 
serendips, each mask was visually inspected by two of the 
authors (BL and RG) independently to search for addi- 
tional serendips using zspec, a publicly available redshift 
measurement program developed by D. Magwick, M. 
Cooper, and N. Konidaris for the DEEP2 survey. In the 
few cases where an object was found by only one of the 
authors or an object was assigned two separate redshifts, 
the slit was "blindly" re-analyzed by a third author (DK) 
and a consensus was reached on the validity and red- 

® See also .http://astro.berkeley.edu/~cooper/deep/spec2d/ 1 



shift of the serendip by all three of the authors. Once a 
serendip was found by eye and confirmed genuine, and if 
spec2d had not detected it on the slit, a manual extrac- 
tion was performed. This process involved re-running the 
spec2d extraction routine on the two-dimensional spec- 
trum with a centroid and FWHM determined by the spa- 
tial location and extent of the serendips as measured in 
the two-dimensional spectrum. This new extraction was 
then inspected and analyzed using in zspec to determine 
if the extraction window was properly centered and the 
aperture was properly matched to the spatial extent of 
the source. In the cases where a non-targeted object was 
detected by eye and spec2d had correctly extracted the 
spectrum, the one-dimensional spectrum was displayed 
with zspec and, if needed, any modifications to the cen- 
troid and FWHM were done iteratively. The redshift 
was determined by guessing the wavelength range of a 
feature (typically 3727 A [Oil], 3968 A CaH, 3934 A 
CaK, 4861 A H/3, 5007 A [OIII], or 6563 A Ha), which 
allows zspec to determine the best-fit redshift through 
an iterative minimization algorithm. All serendips 
found in the C11604 spectral data were found through 
visual inspection, only 30% of which were also detected 
and extracted by spec2d. The small fraction of serendips 
detected by spec2d is not surprising as most galaxies dis- 
covered serendipitously were faint emission-line objects 
and spec2d requires either a continuum or several bright 
emission features to recognize and extract the spectrum 
of a second object on the slit. 

Of the 167 serendips found in this manner, 122 were 
associated with the previously mentioned lower redshift 
(z < 1 for our spectral setup) nebular emission or stel- 
lar absorption lines. The remaining 45 objects were as- 
sociated with either (a) low signal-to-noise ratio (S/N) 
features making a redshift determination uncertain, (b) 
definite features obscured by poor sky reduction or other 
instrumental issues, or (c) a single feature, which in the 
absence of any other spectral indicators makes redshift 
determination difficult, but not impossible (Kirby et al. 
2007). It is the 39 galaxies which comprise category (c) 
that are of interest for this paper. 

2.2. Survey Volume 

The 12 DEIMOS masks observed in the field of the 
C11604 supercluster subtend a total angular area of 3.214 
arcmin^, significantly smaller than the 200 arcmin^ cov- 
ered by the dedicated IMACS Magellan LAE survey 
of M08 and smaller than even the 5.1 arcmin^ covered 
by MS04. However, these surveys have limited vol- 
ume due to their relatively small coverage in the line 
of sight dimension, with spectral ranges comparable to 
that of narrowband imaging surveys (- 100 A). The 
large spectral coverage (6400 A to 9000 A) of the C11604 
DEIMOS data allows for a competitive survey volume. 
The 12 masks sample a volume of 1.70 x 10'* co-moving 
Mpc'^ between z = 4.26 and z = 6.40, slightly smaller 
than other contemporary blind spectroscopic searches 
for LAEs (4.5 x 10^ comoving Mpc^ M08; 6.9 x 10"* 
comoving Mpc'^, S08). However, this volume still does 
not approach the volume covered in narrowband imag- 
ing searches for LAEs such as LALA (7.4 x lO^Mpc^; 
Rhoads et al. 2000; Rhoads & Malhotra 2001), the 
Subaru Deep Field Search (- 1 x lO'^Mpc^; O03), the 
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Subaru XMM - Newton Deep Survey (hereafter SXDF) 
(- 1 X 10^ Mpc^; O08), or the search for LAE galaxies in 
the COSMOS field (~ 1.7 x lO^Mpc^; Murayama et al. 
2007). Though surveying a volume significantly smaller 
than that of narrowband imaging searches, the C11604 
data has the advantage of probing much deeper in the 
luminosity function than such surveys, with a limiting lu- 
minosity of Lii„i — 4.6 X lO"'^ ergs s"""^ at z ~ 5, an order 
of magnitude dimmer than those of narrowband imaging 
surveys (L^m = 4 x 10*^ ergs s~^ , Rhoads & Malhotra 
2001; Lum = 3 X 10''^ ergs s"! , O08; L^m = 6.3 x 10^2 
ergs s""'^ , Murayama et al. 2007). The limiting Ly-a 
luminosity varies slightly (5%-10%) from mask to mask 
due to different integration times and seeing conditions; 
however, the limiting luminosity of Lum = 4.6 x 10^^ ergs 
s~^ represents the brightest limiting luminosity at z '--^ 5 
of all 12 masks, meaning that an LAE at z < 5 with this 
luminosity would be detected in all masks as long as it 
fell relatively close to the center of a DEIMOS slit. 

Three LAEs were detected by M08, representing the 
first successful detection of LAEs by a dedicated spec- 
troscopic survey. Given the survey volume of MOB and 
the range of luminosities found in their survey, it is rea- 
sonable to assume that to detect at least one LAE with 
L> a survey volume of 1.5 x lO'' Mpc'^ is needed. This 
is consistent with the non-detections of T04 and MS04, 
which covered 6.13 x 10^ Mpc^ and 1.1 x 10^ Mpc^, re- 
spectively, and were sensitive to this depth. This limit, 
which excludes the effects of sample variance or any evo- 
lution in the LAE luminosity function between z = 4.26 
and z = 6.4, places our survey right at the volume thresh- 
old necessary to detect a single LAE with L > . 

To calculate the volume of the survey from the entire 
observable redshift range of the DEIMOS masks is, how- 
ever, an overestimate; sky emission features render spec- 
tral regions of the data essentially unusable, necessitating 
bright line fluxes in order to exceed the sky noise. It is 
also tempting at this point to make a correction for the 
angular area of the slit lost by placing a relatively large 
lower redshift object (the targeted galaxy) in the center 
of each slit. However, as discussed in Section 3.2 this 
portion of the slit is not rendered unusable by the target 
galaxy, as we find many serendips and nearly one-third 
of our LAE candidate population at positions coincident 
with the spatial location of the targets. While it is ex- 
tremely likely that the physics governing the observed 
luminosities at these locations differ from serendips dis- 
covered at other positions along the slit (the two most 
likely physical mechanisms are discussed briefly in Sec- 
tion 3.2), this portion of the slit can still be used to 
serendipitously detect galaxies and we therefore include 
it in the calculation of the volume probed by the survey. 
An estimate of the loss due to airglow lines is neces- 
sary, however, and must be done on a slit-to-slit basis 
as the wavelength coverage of each slit is not uniform, 
but depends on the position of the slit along the direc- 
tion parallel to the dispersion on the slitmask, and is 
further compounded by the non-uniformity of the spa- 
tial lengths of the slits. In order to properly account for 
the fractional volume lost by bright sky emission lines, 
we adopt an approach similar to the one taken in SOS. 
For each two-dimensional slit flle, the wavelength value 
of each pixel was determined from the spec2d wavelength 
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Fig. 1. — The top panel shows differential volume (per unit 
redshift) as a function of wavelength for our 12 DEIMOS masks 
in the C11604 field. The vertical lines represent regions straddling 
bright night sky emission features. The lower panel shows the 
cumulative volume of the survey as a function of wavelength as 
corrected for the volume lost by the night sky emission lines (red 
solid line) and the uncorrected volume (blue dashed line). 

solution. Every pixel that was within ± 2a (calculated 
from the FWHM 1200 1 mm^^ resolution) of any bright 
night sky emission line was considered unusable. The 
high resolution of the 1200 1 mm^^ DEIMOS data al- 
lows for minimal losses in usable volume, losing only ± 
1.7 A around each airglow line. Figure [T] shows the us- 
able elements of the data in the spectral dimension as 
well as the cumulative volume covered by the survey as 
a function of increasing wavelength. The volume calcu- 
lated in this manner was 1.365 x 10** co- moving Mpc^, 
~ 20% smaller than that determined by the more naive 
calculation. 

2.3. Flux Calibration 

The DEIMOS spectra were flux calibrated using a flfth- 
order Legendre polynomial flt to time-averaged DEIMOS 
1200 1 mm~^ observations of spectrophotometric stan- 
dard stars^ taken between June 2002 and September 2002 
(see Figure [2]). While the response is known to vary as 
a function of time^, it is a relatively small effect under 
photometric conditions (^5%-10%). As most of our data 
were taken under photometric conditions, we can safely 
ignore this variation. The throughput correction for each 
pixel is: 



Ci Dihc 



r ergs ^ 

cm^ s X' nU9^6x,itexpXc,i, 



(1) 



where Ci are the raw counts in the ith pixel, Di is the 
throughput correction at the central wavelength of the 
ith pixel, 449 is half the effective Keck II mirror aperture 
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Fig. 2. — A fifth-order Legendre polynomial fit to measured val- 
ues of the throughput of DEIMOS for our spectral setup as a func- 
tion of wavelength. This throughput includes loss from other op- 
tical elements and atmospheric transmission. For comparison the 
wavelength range of C11604 LAE candidates found in the data is 
plotted below the throughput curve. The wavelength range of LAE 
candidates encompasses the area of highest instrumental through- 
put. 



in centimeters, S\^i is the plate scale in the ith pixel in 
A pixel" ^, texp is the effective exposure time^, and Ac is 
the central wavelength of each pixel in A. 

The accuracy and precision of the throughput correc- 
tion was checked in the following way. For each high qual- 
ity target galaxy at the redshift of the supercluster, the 
spectrum was multiplied by a fit to the Sloan Digital Sky 
Survey (SDSS) i' filter curve using a quadratic interpola- 
tion to match the wavelength grid of each DEIMOS spec- 
trum. Targets were chosen because they were centered 
widthwise on the slit (serendipitous detections could fall 
anywhere on the slit) and supercluster members were 
chosen because the range of half-light radii was well de- 
termined from the ACS imaging. 

A simulation was run in order to account for losses of 
light due to the finite spatial extent of the slit. Galaxies 
were simulated with exponential disk luminosity profile, 
half-light radii ranging from 0.34" to 0.6", based on val- 
ues measured from ACS F814W data. For each simulated 
galaxy, the light profile was convolved with a Gaussian 
of FWHM comparable to the average seeing conditions 
under which our data were taken (0.9"). A slit of width 
1" and length 6" was then placed on the galaxy, with 
the central part of the galaxy coincident with the cen- 
tral location of the slit. The total flux inside the slit 
was calculated for each simulated galaxy, with the slit 
throughput defined as the ratio of this quantity to the 
total flux in the absence of a slit. This slit throughput is 
plotted as a function of half light radius {r^) and seeing 
in Figure [3] In addition, a similar simulation was run to 
determine the slit throughput as a function of position 
from the slit center under a variety of different seeing 
conditions. Since we are most interested in this effect for 
LAE galaxies, an object with rh = 0.2" was used in the 
simulation, representing a reasonable limit to the sizes of 
large LAEs (see Overzier et al. 2006 or Venemans et al 
2005). 

Despite the functional dependence of slit-loss on the 
object's half light radius, the dependence is not particu- 

^ The effective exposure time is 3600s, as the spectra are nor- 
malized to counts/hour. 
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Fig. 3. — Slit throughput (lo) plotted as a function of half light 
radius (r^) for a variety of different seeing conditions assuming the 
object is placed at the center of the 1" slit (top panel), representing 
an absolute lower limit on the amount of flux that must be lost by 
any galaxy when observed with our spectral setup. The lower panel 
plots the loss of flux as a function of position along the minor axis 
(perpendicular to the spatial axis) of the slit for a galaxy with a 
half light radius of large LAEs (r(j=0.2") for a variety of different 
seeing conditions. While the slit throughput has a moderately weak 
dependence on seeing and half light radius, the dependence on slit 
position is strong, falling off steeply when the object's position is 
more than 0.4" from the central position of the slit. 

larly steep. For objects with rh < 0.4" the dependence 
is essentially linear. Thus, an average slit loss (1 - slit 
throughput) of 0.4 was adopted to correct each spectrum. 
Adopting an average slit loss correction was essential for 
the significant portion of DEIMOS objects which fall out- 
side the coverage of the ACS mosaic and have no reliable 
half light radius measurements. 

The flux density observed in the i' bandpass for each 
spectrum is: 



Sr=0 fx,iS\,iS\.iXc 
v^n Sx.iSx,i 



(2) 



where the sum is over the n DEIMOS pixels that fall 
within the i' bandpass and Sx,i is the i' transmission as 
a function of wavelength. The AB magnitude of each 
spectrum in the i' band was then calculated by: 



« AB,. 



-2.5logih)-A8m--fsec{z), (3) 



with 7 being the airmass term for Mauna Kea^*^. This 
spectral i' spec magnitude was then compared to our Palo- 
mar Large Format Camera (LEG; Simcoe et al. 2000) 
photometry (see G08 for details). Since the slit posi- 
tions were determined from the LEG imaging, there were 
cases where there were noticeable (> 1") positional er- 
rors. Thus, galaxies not centered or absent from the 
slit or those with photometric flags were removed from 
the sample. The derived spectral magnitudes of the re- 
maining galaxies are plotted against the LEG photomet- 
ric magnitudes in Figure [H 

http:/ /www. cfht.hawaii.edu/Instruments/ObservatoryManual/ 
CFHT_ObservatoryManual_(Sec_2).html 
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Fig. 4. — LFC i' magnitudes plotted as a function of magnitudes 
derived from the flux calibrated DEIMOS spectra (top panel) and 
as the difference between the spectral and LFC magnitudes (bot- 
tom panel) for all high quality (Q>3) spectra in the C11604 su- 
percluster. The best-fit relations arc overplotted. The systematic 
offset between the spectral and LFC magnitudes in the brighter 
end of the lower panel is most likely due to an underestimation of 
slit losses or non-photometric considerations. This is not problem- 
atic as the offset drops to zero at the fainter end of the plot, the 
region that our LAE candidates populate. The large rms in the 
lower panel represent real uncertainties in flux calibration of the 
data, a trend that does not improve with decreasing brightness. 

The rms scatter of the spectral magnitudes between 
i' — 19.5 and i' — 25 is 0.49 magnitudes, corresponding 
to an ~ 60% uncertainty in any absolute flux measure- 
ment. While the range of magnitudes are brighter than 
the average magnitude (or magnitude limit) of the LAE 
candidates in our sample, we adopt this rms as being re- 
flective of the uncertainty in Lya line flux measurements. 
In addition, the spectral magnitudes are systematically 
fainter on average by 0.48 magnitudes (lower panel of 
Figured]). While this offset also corresponds to a bias of 
~ 60% for absolute flux measurements, this is less of a 
concern than the rms scatter for several reasons. First, 
the trend in the systematic offset as a function of mag- 
nitude tends towards zero at fainter magnitudes. If a 
magnitude-size relation is assumed for our target galax- 
ies, the observed trend suggests that any offset comes 
from underestimating slit losses for the brighter target 
galaxies. Since we are interested in Lya line fluxes, emis- 
sion which originates from host galaxies that have typi- 
cal i' magnitudes fainter than our dimmest target galaxy 
(i.e., i' > 25.2), this systematic will not adversely affect 
our measurements. While we use a slit throughput of 
0.8 (see the following section) when calculating the line 
fluxes for the purposes of deriving LAE properties such 
as EWs or SFRs, a full slit loss simulation used in calcu- 
lating the luminosity function is undertaken in Section 
5.4. Finally, we approach all measurements from the 
bottom — ,i.e., erring on the side of underestimating the 
true flux of the galaxies so that our measurements will be 
a strict lower limit to compare with other surveys. We 
therefore ignore this systematic and include only the rms 
error when calculating line fluxes. 

2.4. Line Flux Measurements 



For each single emission line galaxy, the one- 
dimensional spectrum was inspected, and three band- 
passes were chosen to measure the emission line flux. 
The first bandpass encompasses the entirety of the emis- 
sion line, avoiding any instrumental or reduction arti- 
facts. The other two bandpasses were chosen to be rela- 
tively sky line free regions blueward and redward of the 
emission line, as close to the emission line in the disper- 
sion dimension as the data would allow, set to a minimal 
width of 20 A. A linear model was fit to each spectrum in 
the blueward and redward bandpasses to mimic the con- 
tinuum throughput. The model parameters were fit with 
a minimization routine, with the associated errors cal- 
culated from the covariance matrix. While a continuum 
fit was typically unnecessary for LAE objects, as the as- 
sociated background was formally consistent with zero 
in most cases, the above procedure was implemented to 
accurately measure the line flux of low-z single-emission 
line galaxies used as a comparison (see Section 3). 

The resulting model background was subtracted from 
each spectrum in the emission line bandpass, with the 
total flux in each bandpass measured by: 



^Lya [ 



ergs 



1=0 



(fA,i(5Ai-B(Ai)) 



1 



(4) 



where B{Xi) is the model at each wavelength, 6Xi is 
the size of the pixel at each wavelength, cosut is the slit 
throughput, and f\ i is defined in Equation [T] The slit 
throughput used in the calculation of the Lya line fluxes 
was set to 0.8, appropriate for a target galaxy with a 
half-hght radius of 0.2" in 0.9" seeing. As most LAE 
candidates are not in the middle of the slit (as a target 
would be) and since the slit-throughput function remains 
below 80% for galaxies centered on the slit for all but the 
smallest half-light radii (r/j < 0.1"), the flux measured 
in this way still represents a lower limit to the true flux 
coming from the galaxy. Tables [T] and [2] list the name, 
redshift (assuming the line is Lya), right ascension and 
declination (assuming the serendip is at the center of 
the slit widthwise), the confidence class, line flux (min- 
imally corrected for fiux losses due to the slit as in the 
above equation), line luminosity, measured or 3<t limiting 
magnitudes, the EW of the Lya line, and the observed 
wavelength of each LAE candidate. 

The associated errors for each fiux measurement were 
derived from a combination of (response corrected) Pois- 
son errors from each spectrum and the errors associated 
with the background model, as well as the flux calibra- 
tion error discussed in the previous section. There can 
also be significant systematic errors associated with the 
bandpass choices. Limiting the size of the emission line 
bandpass can significantly underestimate the true line 
flux, while an overextension of the limits can introduce 
significant noise into the measurement. A select group of 
galaxies, spanning the dynamic range of the spectra mea- 
sured in this manner, were analyzed in order to estimate 
the magnitude of this error. In all cases the system- 
atic errors derived for a "reasonable" range of bandpass 
choices were completely dwarfed by Poisson errors. 

2.5. Flux Limit and Spectral Completeness 

Since our search depended almost entirely on human 
detection of sources, accurately quantifying the com- 
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pleteness limit of the objects detected is more difficult 
than in searches that use automatic peak finding algo- 
rithms. The human eye, while being very good at dis- 
criminating between spurious and real detections and at 
finding irregularities in data (serendips in our case), is 
subject to a variety of effects which are difficult to quan- 
tify. To roughly quantify our completeness limit we sim- 
ulated one hundred slits, each 55 by 8f92 pixels corre- 
sponding to 6.5" by 2700 A at the DEIMOS f 200 1 mm"! 
grating plate scale. These data were first simulated using 
the noise and background properties measured from ac- 
tual DEIMOS two-dimensional spectra in regions where 
features and poor-sky subtraction were absent. These 
feature- free, artifact-free regions were collapsed into one- 
dimensional spectra using the same method used by 
spec2d in extracting one-dimensional spectra of target 
galaxies. Each of the two-dimensional spectra were pop- 
ulated with flux values that mimicked the properties of 
the real two-dimensional spectra, creating in essence one- 
hundred 6.5" "blank-sky" slits. These simulated blank- 
sky slits were populated with objects that varied in both 
intensity and frequency. For each simulated slit, between 
zero and four objects were placed on the slit, character- 
ized by two-dimensional Gaussians with freely varying 
amplitudes, dispersions in both the spatial and spectral 
dimensions, spatial locations, and central wavelengths. 
Noise was also introduced to each Gaussian to properly 
simulate the counting error associated with observing ac- 
tual galaxies. The slits were populated so that a slit had 
zero objects 50% of the time and between one and four 
objects 50% of the time. In addition, the heights and dis- 
persions of the Gaussians were constrained so that the 
objects would have reasonable flux values, i.e. values 
corresponding to an order of magnitude both fainter and 
brighter than the faintest and brightest single-emission 
line object detected in our data. 

Each of the two-dimensional slits was then analyzed 
by one of the authors (BL) in blind observations using 
a fashion similar to that used for the original data. The 
conditions that were present when observing the original 
slits were re-created to the best of our ability (e.g., the 
time spent on each slit, the method of looking for de- 
tections, the software used). For every simulated object 
detected in the two-dimensional slits, a one-dimensional 
spectrum was created using methods similar to spec2d. 
A catalog of generated objects was compared to the cat- 
alog of objects detected by eye and the remaining objects 
that went undetected in the data were then similarly ex- 
tracted. If we set the completeness limit at the faintest 
object detected nearly 100% of the time, this limit cor- 
responds to objects with significances between 105cr and 
111(7 in the two-dimensional data, or a one-dimensional 
significance of 7a. This significance translates to a com- 
pleteness limit of 1.9x 10^^^ ergs s^-^ cm~^ for a 7200s ex- 
posure time, decreasing slightly for our masks with longer 
integration times. This completeness limit is consistent 
with the line flux analysis done in Section 3.3 (see Figure 
[Hand associated discussion), suggesting that this limit is 
close to the actual completeness limit of the survey. 

3. EMISSION LINE TESTS 

The large spectral coverage and moderately high res- 
olution of DEIMOS give us a distinct advantage over 
narrowband imaging searches for LAEs or searches with 



small spectral coverage, as we are able to differentiate 
the Lya line from other emission lines that are typically 
confused for it. The lines which are the most prevalent 
contaminants in searches for Lya emission are the 3727 
A [Oil] doublet, [OIII] at 5007 A, H/3 at 4861 A, or Ha 
at 6563 A. 

The most insidious contaminant in many LAE surveys 
is the [Oil] doublet (rest frame separation 2. 8 A). For our 
spectral setup this line would be observed at a redshift 
of 0.71 < 2 < 1.41 and is usually resolved with the 1200 1 
mm~^ grating. A small fraction of the [Oil] doublets are 
unresolved due to a combination of galactic rotational 
effects and the slit being oriented along the major axis 
of the galaxy. In this case the [Oil] line can still be dis- 
criminated from Lya by the asymmetry of the line. The 
nebular Lya line is typically characterized by its strong 
asymmetry, with suppression of line flux in the blueward 
portion and, in some cases, an extended redward tail. 
A blended (unresolved) [Oil] line in normal star form- 
ing regions (in the absence of an active galactic nucleus 
(AGN)) exhibits asymmetry opposite that of Lya, with 
an extended tail in the blueward portion of the line (Os- 
terbrock 1989; Dawson et al. 2007). Galaxies emitting 
Ha, H/3, or [OIII], in cases of even moderate S/N, can be 
easily distinguished from Lya by other associated spec- 
tral features. The 5007 A [OIII] line is typically seen 
with 4959 A [OIII] and 4861 A H/3 with varying degrees 
of relative intensities (Baldwin et al. 1981). The 6563 
A Ha line can be identified by two accompanying SII 
lines at 6716 A and 6730 A and two [Nil] lines at 6548 A 
and 6583 A, also with varying degrees of relative inten- 
sity. Many spectra originally classified as single-emission 
line objects were recognized as low redshift interlopers 
through the identification of faint associated lines. 

For the remaining 39 objects that were classified as 
genuine single-emission line objects, several tests were 
performed to further remove any low redshift interlop- 
ers. The Lya line is characterized by a large 1.3-4.5 
mag continuum break blueward of the line due to atten- 
uation of Lya photons by intervening neutral hydrogen 
(H04). Initially, the spectral data were inspected, and 10 
single-emission line serendip exhibiting appreciable con- 
tinuum blueward of the emission feature relative to any 
redward continuum was eliminated as a potential LAE 
candidates. The imaging data were also useful in dis- 
criminating single-emission line serendips in this regard, 
as the photometric filter setup would also, in many cases, 
probe the continuum break across the Lya line (see Fig- 
ure[5]). Each single-emission line serendip detected in one 
or more of the photometric bands was required to exhibit 
a continuum break over filters blueward and redward of 
the line. Since most of these objects are extremely faint 
in the imaging (if they are detected at all), requiring a 
strong continuum break over the emission line is, in al- 
most all cases, similar to requiring that the object drop 
out of any band blueward of the emission line. Our bluest 
LFC and ACS bands, r' and F606W, are situated so that 
either would pick up a significant amount of continuum 
fiux from any LAE at the bluer end of our detection limit 
(A < 7000 A, 

ZLya ^ 4.75). For objects such as this we 
had to rely on Subaru Suprime-cam V-band data to dis- 
criminate between potential LAEs and low-z interlopers. 
Any galaxy detected in the V-band data was excluded 
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Fig. 5. — Synthetic Lya spectra at three different redshifts 
spanning the range of redshifts of our sample overlaid on the four 
ground-based filter transmission curves (top three panels) and two 
ACS filter transmission curves (bottom three panels) for which 
we have coverage. The continuum break over the Lya line is mod- 
eled based on narrowband imaging measurements of the continuum 
break {H04). For lower redshift LAEs, the only band able to detect 
the continuum break is the ground-based V band. At the redshift 
of most of our sample (z ~ 4.8), the Lya line just passes the cover- 
age of the ground-based r' band and lies just at the red end of the 
F606W band. At higher redshifts the Lya line is comfortably situ- 
ated in the F814W band and at the red end of the ground-based i' 
band, giving us significantly more power to discriminate between 
interlopers and genuine LAEs at these redshifts. The throughput 
of the V band is scaled down and the throughput of the z' band is 
scaled up for clarity. 



as an LAE candidate due to the relative shallowness of 
the image (see Section 4.1 for details on the depth of 
the photometry). All but two of 22 single-emission line 
galaxies that were eliminated as potential LAE candi- 
dates through the above tests failed the continuum break 
test. The two single-emission line low-z interlopers that 
did not fail this test were among the 10 galaxies that 
failed the spectral continuum break test. In addition, 
each single-emission line serendip that was detected in 
the photometric data was also inspected visually, and 
any objects with large (> 2") angular extents were clas- 
sified as low-z interlopers. Six of the 22 single-emission 
line \ow-z interlopers were rejected by this test, although 
in all cases these galaxies had failed at least one of the 
two previous tests. 

Of the original 39 single-emission line cases, 17 objects 
survived the previous tests. A small subset of these ob- 
jects were insensitive to these tests, as the single-emission 
line object was superimposed spatially in the spectral 
data with either the target or another serendip. In such 
cases, the single-emission lines were checked against a va- 
riety of nebular emission lines at the redshift of the super- 
imposed target or serendip to verify that it could not sim- 



ply be an unusual emission feature coming from the same 
galaxy. In many cases, however, it was clear from the 
morphology or positions of the lines that the two emission 
features originated from two separate sources. In some 
cases the two superimposed objects were resolved in the 
ACS data, and the continuum break test was used on one 
or both of the galaxies, depending on whether the iden- 
tity of the single-emission line source was certain. More 
frequently, however, the two objects remained unresolved 
in the HST data, so we include them in our sample. Of 
the 17 objects that survived the original single-emission 
line tests, all 17 passed the tests described above. These 
17 objects comprise our LAE sample (see Figures [T91I2T1) . 

3.1. Line asymmetry 

Another discriminator used on the individual single- 
emission line spectra was a computation of the wave- 
length asymmetry parameter (Dawson et al. 2007). 
Briefly, the asymmetry parameter, aA, is defined as: 



1 



A. -A 



10, f) 



Al 



A. 



(5) 



where Ac is the central wavelength of the emission, de- 
fined as the point of maximal flux in the line profile, and 
Aio,r and Aio,& are the wavelengths where the flux first 
exceeds 10% of the peak flux redward and blueward of 
the line, respectively. This diagnostic can be used to fur- 
ther discriminate single-emission lines that exhibit stan- 
dard Gaussian (Voigt) profiles such as II/3, Ha, [OIII] 
(1/aA ~ 1), or a blended 3727 A [Oil] doublet [l/ax > 
1) from a higher redshift Lya line that exhibits strong 
asymmetry in the opposite direction. While this test can 
be a useful diagnostic in a statistical sense, an asym- 
metry parameter of I/oa > 1 was not a strong enough 
constraint to rule out an object as an LAE candidate 
if it had passed all the previous tests. This is because 
several processes (instrumental broadening, local under- 
densities of HI regions, etc.) can cause the LAE emission 
to appear symmetric. Conversely, an object which had 
failed one or more of the above tests was not reclassified 
as a potential LAE based on an unusually high (I/aa < 
1) asymmetry parameter, as low redshift lines can, under 
rare circumstances, exhibit strong redward-skewed asym- 
metry (see for example object D21 in MS04). Therefore, 
this diagnostic was used only to discriminate between 
high quality LAE candidates and poorer quality can- 
didates, rather than distinguishing genuine LAEs from 
interlopers. Figure [6] shows a histogram of the inverse 
of the asymmetry parameter of the known lower redshift 
single-emission line objects, a population of blended [OH] 
emitters (confirmed by other associated lines present in 
the spectrum), and our 17 LAE candidates. The objects 
clearly separate out; the LAE candidates primarily oc- 
cupy the high asymmetry (low inverse asymmetry) por- 
tion of phase space, the low-z interlopers are distributed 
around unity (symmetric) , and the [OH] galaxies are pri- 
marily situated in the region of phase space opposite that 
of the LAE candidates. In fact, all but three LAE candi- 
dates (all Quality 1; see below) have inverse asymmetry 
parameters less than 0.75. 

3.2. LAE Confidence Classes 
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Fig. 6. — The inverse of the asymmetry parameter measured for 
the emission lines of three different populations: known blended 
[Oil] emitters (0.7 < z < 1.4), low-z single-emission line galaxies, 
and our 17 LAE candidates. This parameter easily discriminates 
even blended [Oil] emission lines, which are strongly skewed toward 
high values of the inverse of the asymmetry parameter, from LAE 
lines which typically exhibit values much less than unity. The 
low-z single-emission line galaxies are less easily discriminated by 
this test, with a distribution that is centered around symmetric 
(l/a;^=l) line profiles and wings that extend into the phase space 
of the other two populations. All LAE candidates with higher 
inverse asymmetries {1/ax > 0.75), including one very high value, 
are low-quality candidates. 

Each of the 17 LAE candidates was assigned a quahty 
class. Quality classes are assigned to LAE candidates 
in a fashion nearly identical to that of SOS and are de- 
fined as follows: Quality 1 objects pass all of the above 
tests, but show no additional indicators of being gen- 
uine LAEs. Objects which are Quality 1 do not exhibit 
any asymmetry (or exhibiting blueward-skewed asym- 
metry) in their line profiles and are non-detections in 
all photometric bands. These objects are our least se- 
cure candidates, nearly equally likely to be low lumi- 
nosity foreground galaxies as LAEs. Quality 2 and 3 
objects all similarly pass the interloper tests but also 
show strong asymmetric line profiles. A few of these 
objects are detected in one or more photometric bands, 
further increasing our confidence in these objects as gen- 
uine LAEs, but it is the asymmetric line profile which 
is the defining characteristic of the higher confidence 
classes. Both Quality 2 and Quality 3 candidates rep- 
resent our highest level of confidence that an object is 
a genuine LAE. However, Quality 2 objects are super- 
imposed with a target or another serendip spatially on 
the slit. Thus, the flux measurements of the Quality 2 
objects could be signiflcantly dimmed by extinction from 
the interstellar medium (ISM) of the foreground galaxy 
or boosted through galaxy-galaxy lensing. This addi- 
tional, unknown component of the uncertainty makes it 
necessary to exclude Quality 2 galaxies from certain parts 
of the analysis. 

3.3. Flux and Redshift Tests 

The tests discussed in the beginning of Section 3 can 
only be used to rule out objects as LAEs, not to prove 
that any particular object is definitively an LAE. The 



tests in the following two sections explore the statisti- 
cal similarities or differences between LAE candidates 
and the low-z interlopers, giving us further confidence 
that the LAE candidates represent a unique and sepa- 
rate population. 

3.3.1. Effective Redshift Test 

First we compare the observed wavelengths of the 
single-emission lines in the low-z interloper population 
to the observed wavelengths of the Lya lines in the LAE 
candidates. The low-z single-emission line interlopers are 
comprised of some combination of [Oil], H/3, [OIII], and 
Ha emitters and therefore cannot be given definite red- 
shifts. Following the analysis done in SOS, we have recast 
the low-z interlopers in terms of an effective redshift: the 
redshift that the object would have if the line were Lya, 

such that Zeff,Lya = (Aem/1215.7 - 1). 

The idea of this test is that the low-z single-emission 
line interlopers, if they truly are comprised of a mix of 
the aforementioned lines, should be, in the absence of 
any instrumental effects, equally distributed in effective 
redshift (wavelength) space. An object at a redshift of 
z = 0.35 emits the 5007 A[OIII] line at Xobs = 6759 
A and 6563 A Ha at Xobs = 8860 A, both of which 
could mimic single emission lines under a variety of dif- 
ferent conditions. These effects could be (1) instrumen- 
tal: the placement of the slit on the slit mask truncating 
either the blue or red end of the CCD response; (2) at- 
mospheric: a bright night sky line masking the second 
emission line; or (3) a result of galactic processes: a low 
level AGN which exhibits strong [OIII] emission but lit- 
tle to no Balmer emission, or a starburst galaxy having 
strongly suppressed forbidden transitions relative to the 
strength of the Balmer lines. In any of these cases, the 
chance is more or less equal that the single-emission line 
galaxy will show up as the blue or the red emission line. 
The redshift distribution of the LAE population should 
be strongly biased towards the lowest redshifts to which 
we are sensitive, as we probe successively shallower in 
the luminosity function as the LAEs move to higher red- 
shifts. Thus, if the LAE population represents a truly 
different population than the low-z single-emission line 
interlopers, the redshift histograms should differ signifi- 
cantly. 

Figure [7] shows the comparison in effective redshift 
space between the 22 low-z interlopers, the 4 Q=l and 
the 13 Q=2,3 LAE candidates. The low-z single-emission 
line interlopers are more or less evenly distributed across 
Zeff,Lya with two important exceptions. There are no 
interlopers shortward of Zeff,Lya = 4.4, possibly due to 
the prevalence of Ha as the unknown single emission line 
in the interloper population. The rest wavelength of the 
Ha line has Zeff,Lya = 4.398 so if the interloper popula- 
tion does consist primarily of Ha emitters, few galaxies 
would be seen blueward of this limit. Another reason 
for this drop-off in detections could be the significant 
drop in DEIMOS sensitivity blueward of ^ 6600 A for 
our spectral setup. The second drop-off in detections oc- 
curs at Zeff^Lya > 6.1, most likely due to the significant 
decrease in DEIMOS sensitivity and the decrease in sig- 
nificant sky line free spectral windows redward of 8700 
A (see Figure [2]). 

The LAE population is strongly peaked towards the 
low end of our redshift sensitivity. A very noticeable 
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Fig. 7. — Redshift histogram of the high quality (Q=2,3) and 
low quality (Q=l) LAE candidates. A strong peak can be seen at 
2 ~ 4.8 as well as less pronounced peaks at 2 ~ 4.4 and z ~ 5.7. 
Conversely, the overplotted low-z interloper population (plotted in 
terms of z^ff i^y^; see Section 3.3.1) is distributed nearly symmet- 
rically around z ~ 5.1, with a much slower falloff at high redshifts. 

peak exists at z ^ 4.8, which may represent a real clus- 
tering of the LAE population in projection space or could 
simply be an artifact of the sensitivity issues discussed 
in the previous paragraph, as the DEIMOS sensitivity 
peaks at ~ 7000 A for our setup. More likely, it is some 
combination of these two effects (see Section 5.2 for a 
discussion). The Q=l LAE candidates, which are our 
least secure candidates, are surprisingly consistent with 
our higher confidence Q=2,3 population, also peaking 
around ZLya ~ 4.8. There are two Q=2,3 candidates at 
z ~ 5.7 which are unexpected, given our prediction that 
the LAE population should be strongly peaked towards 
the low-redshift end. 

3.3.2. Lme Flux Test 

The second of these tests explores the possibility that 
the LAE candidate population represents a lower lumi- 
nosity subset of the single-emission line interlopers. A 
majority of single-emission line interlopers were ruled 
out by broadband detections, i.e., not exhibiting a suf- 
ficiently strong continuum break over the feature to be 
plausibly identified as Lya. All of the single-emission 
line interlopers were detected in the photometry. Con- 
versely, the majority of the LAE candidate population 
were not detected in any of the three LFC bands nor the 
two ACS bands. Thus, the LAE candidate population 
clearly represents a class of objects that are significantly 
dimmer in continuum luminosity. If the LAEs are truly 
drawn from the same population as the low-z interlop- 
ers, their line luminosities should similarly scale down. 
This test provides a quantitative statistical tool to dif- 
ferentiate the LAE candidates from the lower luminosity 
tail of the single-emission line interloper luminosity func- 
tion. This test is not sensitive to the case where the LAE 
candidates represent a population of dwarf starbursting 
galaxies with higher line luminosity relative to their con- 
tinuum brightness (Fricke et al. 2001; Guseva et al. 2003; 
Kehrig et al. 2004; Izotov et al. 2006). 

Figure [8] shows the comparison between the line fluxes 
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Fig. 8. — Histogram of line fluxes uncorrected for slit losses 
of our LAE candidates and confirmed low-z single-emission line 
objects. While the two distributions overlap, the low-2 interloper 
population is characteristically dimmer than both the high {Q=2,3) 
and low {Q=l) quality LAE candidates. The difference in the 
observed distribution suggests that the LAE candidates are not 
primarily comprised of low luminosity single-emission line objects 
at low (z < 1.4) redshifts. The vertical long dashed line represents 
our adopted completeness limit of 1.9 X 10~^* ergs cm~-^ 

of the single-emission line interlopers relative to the LAE 
candidates. The Q=2,3 LAE candidates are, on aver- 
age, brighter than the single-emission line interloper ob- 
jects, with the mean line flux about 0.5 dex higher than 
the interloper population. The average magnitude of 
the interloper population in the band which best sam- 
ples the continuum emission near the emission feature 
is 23.5 mags. In contrast we can adopt the LFC i' 3a 
limit of 24.3 mags as the upper bound on the continuum 
flux of LAEs that are not detected in the photometry 
(a conservative limit as many of the candidates are un- 
detected in the ACS images which have a 3a depth of 
~ 26 mags). This limit on the continuum flux requires 
the LAE candidates, if they are instead low-luminosity, 
low-z interlopers, to have line equivalent widths (EWs) 
at least 10 times greater than the average EW of the 
known single-emission line interlopers (5.4 A). Such high 
EWs are certainly plausible in dwarf galaxies undergoing 
a starbursting event where the EWs of Ha (usually the 
strongest lines in optical starbursting spectra) are the 
range 50-150 A (Kennicutt 1998, Petrosian et al. 2002) 
and have been measured as high as ~ 1500 A (Kniazev et 
al. 2004, Reverte et al. 2007). However, such objects are 
uncommon, and we would expect to observe other associ- 
ated lines (e.g., [Nil], H/3, [OIII]) in the data, which we do 
not. It is interesting to note that if we adopt a standard 
ratio for log[i([iV//] 6585 k)/L{Ha)] and log[L([0///] 
5007 k)/L{HI3)] of -0.45 for star-forming galaxies (Bald- 
win et al. 1981; Brinchmann et al. 2004; Shapley et al. 
2005; Yan et al. 2006), the bulk of our LAE candidates 
(^60%) are sufficiently brighter than the completeness 
limit so that the associated lines would be detected if 
the emission were instead Ha or H/3. 

The Q=l LAE candidates are essentially identical to 
the fluxes of the low-z interloper population, with a mean 
flux of 4.8x 10^^® ergs s^^ cm~^ as compared to the mean 
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flux of the interloper population of 3.9 x 10^^ ergs 
cni~^. While this similarity may be an indication that 
the Q=l LAE candidates contain at least some \ow-z 
interlopers mixed in with genuine LAEs, it also may be 
misleading. The average upper limit on the magnitude of 
the Q=l candidates in the filter sampling the continuum 
surrounding the emission feature is 24.9, nearly 1.5 mag- 
nitudes dimmer than for the single-emission line counter- 
parts. While the line fluxes of these two populations are 
similar, the EW of the Q—1 LAE candidates would still 
necessarily have to be a factor of 4 higher than the in- 
terloper population. In addition, the Q=l line fluxes fall 
near the completeness limit of 1.9 x 10~^^ ergs s~^ cm~^ 
and near the low-flux tail of the line flux measurements 
of the high quality (Q=2,3) LAE candidates. We would 
expect, independent of the redshift range, an inverse re- 
lationship between the number of detections and the line 
flux down to the completeness limit and a steep falloff in 
detections thereafter. If the Q=l LAE candidates con- 
stitute real detections of genuine LAEs, this would be 
the behavior we observe in the data. Thus, it may be 
that these low quality candidates simply represent the 
fainter flux end of the LAE population, and their lower 
S /N prevents them from reliably being classified as high 
quality candidates. 

3.4. Composite Spectra 

Previous tests focused on measurements of individual 
spectra of galaxy signatures at or near the flux limit, 
making these measurements susceptible to noise effects. 
While we compare the ensemble properties of the galaxy 
populations, which is less sensitive to noise variations 
in the data than the comparison of individual measure- 
ments, an alternative is co-adding of the spectra in order 
to increase the S/N. 

To properly retain the overall spectral properties of the 
constituent objects (e.g., line shapes, resolution, velocity 
dispersions, etc.) and to avoid averaging out faint fea- 
tures, it is necessary when coadding galaxies to determine 
the redshift as accurately or in as consistent a manner as 
possible. While we were able to determine redshifts for 
our interloper sample through the centroiding provided 
by spec2d, the LAE population was problematic because 
of the uncertainty in determining the true peak of the 
line. In the absence of any other knowledge about the 
true profile of the Lya emission in each galaxy (other 
than its asymmetry), the assumption was made that the 
wavelength associated with the peak fiux in each emission 
line profile represented the central wavelength for that 
emission. Lack of knowledge about the shape of the true 
line profile introduces a significant {5z ~ 0.002) absolute 
error in the redshift measurements. However, since this 
measurement is made in a consistent way for each spec- 
trum, the relative error in the redshifts (the important 
quantity for coadding purposes) between any two spectra 
is quite small. Thus, any asymmetry in the original line 
profiles should be preserved through this process. 

Each galaxy spectrum was then "de-redshifted" to its 
rest frame, or, in the case of the single-emission line in- 
terlopers, the effective rest frame (see Section 3.3). Each 
rest frame spectrum was interpolated onto a pixel grid 
of common size, chosen to subsample the lowest (rest- 
frame) pixel scale. The resulting spectra were then added 
together in the following two ways: (1) each spectrum 




Fig. 9. — Composite spectrum of the LAE candidates and the 
known low-z interlopers smoothed with a Gaussian smoothing ker- 
nel of o'=1.5 pixels. Each co-addition is done using luminosity 
weighting (see the text). The dashed lines show Gaussian fits to 
the line profiles. Both the high (Q=2,3) and low (Q=l) quality 
LAE candidates are poorly fit by a Gaussian model. 

was normalized by the galaxy's total spectral fiux (uni- 
form weighting), or (2) galaxies were added together with 
no normalization (luminosity weighting). In both cases, 
the fiux of each pixel in the co-added spectrum was pop- 
ulated using a Poissonian variance weighted mean of the 
pixel values at each wavelength in the individual spectra. 

Figure [9] shows the luminosity- weighted coadded spec- 
trum for three different populations: the high quality 
(Q=2,3) LAE candidates, the low quality (Q=l) LAE 
candidates, and the known low- 2: interlopers. The coad- 
ded spectrum of each set of galaxies was fit with a Gaus- 
sian, with the goodness of fit parameterized by the re- 
duced • As expected, the Gaussian model does a poor 
job at reproducing the observed line profile for the high 
quality LAE candidates {x^ /v = 2.38). Conversely, both 
the low quality LAE candidates {x^ jv — 0.67) and the 
low-z interloper {x^/v = 0.88) population are statisti- 
cally well fit by the Gaussian profile. Despite the statis- 
tical significance of the fit, visual inspection of the low-z 
interloper population shows the line profile to be clearly 
more symmetric than the low quality LAE candidates, 
as should be the case if at least some of the low qual- 
ity LAE candidates are real. Additionally, the best-fit 
Gaussian to the low quality candidates has a FWHM 
of 0.68 A, nearly twice as large as the best-fit FWHM 
of the known low-z interlopers, further suggestive that 
the lower quality LAE population contains at least some 
genuine LAEs. The inverse of the asymmetry parame- 
ter (Section 3.1) was also calculated for the luminosity- 
weighted CO- added spectra of each of the galaxy subsets, 
with values of 0.35, 0.57, and 1.14 for the high-quality 
LAE candidates, low quality LAE candidates, and low-z 
interloper population, respectively. Both of these results 
reinforce the conclusions reached from analyzing the in- 
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dividual spectra: that the high quahty LAE candidates 
probably represent a real population of LAEs while the 
lower quality candidates probably represent some com- 
bination of genuine high redshift LAE galaxies and low- 
z interlopers. The results of these calculations did not 
change significantly if we instead use uniform weighting. 

4. PROPERTIES OF THE CL1604 LYq EMITTERS 

4.1. Photometric Limits 

The broadband photometry associated with the C11604 
data set was designed almost exclusively to select spec- 
troscopic targets for the supercluster at (z) — 0.9, sam- 
pling down to 3tT limits of 24.8, 24.3, and 23.6 in r', i', 
and z' respectively. These magnitudes were calculated 
by measuring the magnitude of a circular object with a 
1" diameter, where each pixel has signal equal to three 
times the sky rms (effectively a circular top hat profile). 
An aperture of 1" was chosen to match the average seeing 
conditions from on Paloniar mountain during our obser- 
vations. 

The depth of these observations are only sufficient to 
probe the continuum luminosities of the most massive 
galaxies at high redshift (z > 4.4). Indeed, only one 
of our LAE candidates (16XR1.72, an object that was 
subsequently picked as a spectroscopic target) was de- 
tected to the depth of these images. The accompany- 
ing archival Suprime-cam observations have a 3cr limiting 
magnitude ¥^^24.0 for the same choice of aperture as the 
other ground-based images. The exact value of this limit 
is unknown due to imperfect photometric calibration, 
though it is probably accurate to ~0.2 mags based on 
comparisons between the measured Subaru magnitudes 
and overlapping fields with precise photometric calibra- 
tion. 

The ACS observations are significantly deeper, reach- 
ing 3cr limits of 26.1 and 25.5 in F606W and F814W 
in most of the pointings and 26.8 and 26.3 in two deeper 
pointings centered on clusters A and B. Photometric lim- 
its in the ACS pointings are calculated for a 0.3" circu- 
lar aperture using the same method as the ground based 
limiting magnitudes. A smaller aperture was chosen be- 
cause of the significant increase in resolution ACS pro- 
vides relative to the ground based images. These 3(t lim- 
its are conservative limits on the depth of our images as 
the differential number counts do not turnover (hereafter 
"turnover magnitude") until magnitudes that are 0.1-0.2 
fainter than the 3cr limits of the ground-based data and 
0.5-1 mags fainter than the limits of the ACS data. Even 
though the ACS data does not overlap the entirety of 
our spectral coverage, only two of our 17 LAE candidates 
(SC2NM1.45 and SC2NM2.61) feU outside the ACS area. 
Despite this, only three of the 15 LAE candidates that 
were covered by ACS pointings were detected in the ACS 
imaging. 

In order to place limits on the broadband photometry 
in the absence of detections, local versions of 3(7 limiting 
magnitudes were measured for each LAE candidate from 
the data using a method similar to the measurement of 
the 3(7 limiting magnitudes for each image. However, 
rather than measuring the rms over a large portion of 
the image, the rms was instead measured in a statisti- 
cally significant region cither at the central location of 
the galaxy (inferred from the spectroscopy, assuming the 
object was at the center of the slit) or near the target lo- 



cation if the object was superimposed spatially with the 
target. For 3cr limiting magnitudes in the ACS images, 
this rms value per pixel (corrected for correlated noise 
from pixel subsampling) was multiplied by the number 
of pixels covered by an object with a circular aperture 
of radius 0.21". This number was motivated by the half 
light radius of LBGs (Steidel et al. 1996b; Ferguson et 
al. 2004) and intentionally designed to overestimate the 
limiting magnitude of such objects; all of the LAE can- 
didates detected in the ACS imaging had detected mag- 
nitudes significantly dimmer than the corresponding 3cr 
limiting magnitude. In addition, the 3fT limiting mag- 
nitudes were measured with the Palomar LFC imaging 
using similar techniques. As before, a circular aperture 
of 1" was used in the LFC calculation, as a typical LAE 
would not be appreciably different spatial extent than a 
point source in the LFC images. Table[T]gives the 3ct lim- 
iting magnitudes of all the non-detected LAE candidates 
as derived from both the ACS imaging (when available) 
and the Palomar LFC imaging. 

4.2. Equivalent Widths 

The EW is typically calculated for the Lya line in the 
following way (Dawson et al. 2004): 

EWiLya) = (6) 
/a(1 + z), 

where FLya is the total line flux in the Lya line and fx 
is the flux density redward of the Lya emission, a for- 
malization that is convenient for measurements of LAEs 
in narrowband imaging surveys. 

Without proper detections of the continuum luminos- 
ity of a majority of our LAE candidates, calculating the 
EW of the Lya line, something that is strongly depen- 
dent on the continuum luminosity, is not possible. In- 
stead we calculate a lower bound on this quantity. For- 
mally, our 3(7 limiting magnitude represents a strict up- 
per (brighter) bound on the continuum flux density. The 
uncertainty in the flux loss in the Lya line due to the 
slit works in the same direction; the total line flux, mini- 
mally corrected for slit losses (see Section 2.4), represents 
a strict lower (dimmer) bound on the line flux. Thus, any 
calculation based on these numbers will represent a very 
conservative lower bound to the EW of the Lya line in 
these galaxies. 

In order for the EW measurement, or a lower bound 
to this measurement, to characterize the intrinsic prop- 
erties of high-redshift LAEs, it is necessary to make some 
correction for attenuation from the intergalactic medium 
(IGM). This attenuation occurs primarily due to reso- 
nant scattering of redshifted Lya photons in intervening 
clouds of neutral hydrogen. As such, only Lya photons 
emitted by galactic components blueshifted with respect 
to the bulk velocity of the galaxy will be affected by 
this dampening. Although there can be, in principle, 
some contribution to the attenuation from intervening 
Helium and metal systems, such contributions are typi- 
cally small in comparison (Madau 1995). The attenua- 
tion to the blueward flux solely from intervening HI re- 
gions was characterized most recently by Meiksin (2006), 
where the fraction of attenuated Lya photons blueward 
of 1215.7 A was given as: 
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fau.Lya,b ^l-exp (-5.8 X 10-4 (1 + z)4-5) (z > 4), 

where the argument of the exponential is the mean Gunn- 
Peterson optical depth for an object at a given redshift. 
Assuming the LAE is rotationally supported such that 
there is no skew in the velocity components of the Lya 
emitting HI regions, the true flux of the Lya line in LAEs 
(for z > 4) is given by: 



F, 



F, 



corr,Lya. 



Lya 



0.5 + 0.5(1 



(8) 



an expression that ignores any dust extinction of the Ly- 
man continuum. While Equation [7] is derived from an 
average of different lines of sight from observed data, we 
use it here to correct on a galaxy by galaxy basis. Though 
making this correction may introduce significant bias to 
the EW measurement of a single galaxy, correcting our 
entire sample produces a distribution which more accu- 
rately reflects the true contribution of star-forming pro- 
cesses in these galaxies. After correcting each galaxy's 
line flux using Equation [H the upper bound of the con- 
tinuum flux density, f\, was estimated. For the bulk of 
our sample which went undetected in the photometry, 
the flux density was estimated with both the Scr limiting 
magnitude in the band encompassing the Lya emission 
and in a band just redward of the Lya emission. For 
the higher redshift galaxies (z > 5.5), we had no bands 
completely redward of the Lya line with sufficient depth 
to make a meaningful estimate of the EW with the LEG 
data (see Figure [5]), as our z' imaging was shallower than 
our other bands. Both cases the LAEs fall within the 
AGS imaging and we therefore use only the 3cr magni- 
tude limit in the F814W filter to estimate the flux density 
redward of the Lya line. 

Since these galaxies are undetected in the AGS or LEG 
data, both the EW measurement from a band encom- 
passing the Lya line and from a band longward of the 
Lya line represent a reasonable approximation to the 
lower limit on the true EW. Although we formally cal- 
culate EWs from the 3cr limits in bands containing the 
line, the true lower bounds to the EWs are characterized 
solely by the EWs calculated from the Scr limiting mag- 
nitudes of bands redward of the Lya line. In addition, 
we have calculated the EWs using the turnover magni- 
tude (see Section 4.1 for deflnition) in the AGS imaging. 
For the shallower AGS pointings this magnitude corre- 
sponds to 26.99 and 26.68 in F606W and F814W respec- 
tively, increasing to fainter magnitudes (27.76 and 28.01 
in F606W and F814W) for the deeper AGS pointings. 
These turnover magnitudes are not to be confused with 
the completeness limits of the images, which must be 
constrained through simulations; however, the turnover 
magnitudes are likely a good approximation to the limit 
at which we are complete for LAE-size objects. While 
the EWs calculated using the turnover magnitudes are 
not strictly lower limits, they serve as more reasonable 
(less conservative) estimates of the lower bound of the 
Lya EW (see Tabled]). 

For galaxies detected in the imaging (either AGS or 
LEG), the calculation was much more direct. While any 
measurement of the EW still represents a lower limit (due 
to the unknown amount of flux loss of the Lya line from 
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Fig. 10. — The distribution of lower limits to the EW of the 
LAE candidates in our sample. Lower limits were calculated us- 
ing line fluxes minimally corrected for slit losses and either the 
detected magnitude or the magnitude of the completeness limit in 
a filter completely redward of the Lyo emission (for z < 4.9) or 
encompassing the Lyo emission (for z > 4.9). 

the slit), the redward continuum flux densities could be 
calculated in a straightforward way from the observed 
magnitudes. These measurements were done, as was the 
case for the undetected objects, for both the band en- 
compassing the Lya emission and a band redward of the 
line emission. The lower limits on these EWs are shown 
in Table [21 quoted at the 95% confidence level. 

Figure [To] shows the distribution of our EWs calculated 
from either the broadband magnitude completeness limit 
or the broadband magnitude of the detection redward of 
the Lya line. The distribution is strongly skewed towards 
very modest values (~25 A) of the EW as compared to 
measurements from other surveys (Dawson et al. 2004, 
II04, O08). This result is not surprising, due to the bulk 
of our sample populating the faint end of the Lya lumi- 
nosity function and to the manner in which we estimate 
the continuum luminosity of candidates undetected in 
the imaging data. Since the continuum luminosity is es- 
sentially independent of Lya luminosity and since the es- 
timated continuum luminosity is most likely significantly 
brighter than the true continuum luminosity of our can- 
didates, the observed distribution may be more reflective 
of the way in which the EW was calculated rather than 
of any intrinsic properties of the LAE candidates. Still, 
the observed EW distribution is comparable to the re- 
suhs of the GLARE survey (Stanway et al. 2004, 2007), 
a comparably dim sample of LAEs, suggesting that there 
may be inherent properties of low- luminosity LAEs which 
contribute to our observed distribution. 

4.3. Star- Formation Rates and Star- Formation Rate 

Density 

For each LAE candidate the star formation rate (SER) 
was calculated by: 

SFR [M^yr-^] = 9.5 x IQ-^^ L^orriLya), (9) 

where Lcorr{Lya) is the line luminosity based on the 
calculation in Section 2.4 and Equation 8. The constant 
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of proportionality in Equation [9] is derived from the Ha 
relation used in Kennicutt (1998) which assumes contin- 
uous star formation with a Salpeter IMF and the ratio of 
Lya to Ha photons calculated for Case B (high [T(Lya) 
^ 10"*] optical depth) recombination (Brocklehurst 1971) 
for an electron temperature of Tg = lOOOOK and solar 
abundance. This corrected SFR still represents a lower 
limit to the intrinsic SFR of the galaxy due to the un- 
known amount of slit loss and due to the fact that we do 
not correct for dust extinction. Tables [T] and [2] list the 
calculated SFRs for all 17 LAE candidates. The typical 
LAE candidate galaxy in our sample forms stars at a rate 
of 2-5 M0 yr"'^, on the low end of SFRs found in other 
samples. 

As a consistency check, we also calculated SFRs from 
the rest-frame ultraviolet continuum for the three LAE 
candidates that had detectable continuum in bands red- 
ward of the Lya line. In no cases was this measure- 
ment possible from the spectra as the continuum strength 
in the spectra redward of the Lya line was not strong 
enough for a reliable measurement. Each of the UV SFRs 
were derived from the flux density of the three photo- 
metrically detected LAEs, measured from their F814W 
magnitudes and converted using the Madau et al. 1998 
formula: 

SFRuv [Moyr^i] = L,,uv/{S x 10^^), (10) 

where J^u,uv is a luminosity density calculated at ^^1500 
A. Since the three LAEs detected in the photometric 
data are at different redshifts, the effective rest-frame 
wavelength for the F814W filter changes slightly. How- 
ever, we made no attempt to K-correct the observed flux 
densities as the effective rest-frame wavelength is always 
within 80 A of 1500 A and because the spectrum of LAEs 
are relatively flat in the UV. Table[2]lists the UV-derived 
SFR as well as the absolute UV magnitude for each of 
the three LAE candidates detected in the imaging data. 
The values of the SFR derived in this manner are very 
similar to the SFR derived from the Lya line, with the 
exception of one of the candidates (FG1.20), suggesting 
that slit losses for two of these candidates is minimal 
(not surprising since one of the objects, 16XR1.72, was 
targeted) . 

Another quantity of interest for any population of 
high redshift galaxies is the star formation rate den- 
sity (SFRD), which can be used to determine the on- 
set of reionization in the universe. Observations of very 
high redshift (z > 6) LAEs and LBGs (Malhotra & 
Rhoads 2004, hereafter MR04; Kashikawa et al. 2006, 
hereafter KG6; Shimasaku et al. 2005; Taniguchi et al. 
2005; Bouwens et al. 2003, 2004b, 2006; Bunker et 
al. 2004), Gunn-Peterson troughs in very high-redshift 
quasars (Becker et al. 2001; Djorgovski et al. 2001; 
Fan et al. 2006), and optical depth measurements from 
the Wilkinson Microwave Anisotropy Probe (WMAP; 
Spcrgel et al. 2007; Hinshaw et al. 2009) all suggest that 
the universe was reionized many 100s of Myrs prior to the 
observational epoch of our sample. The dominant pop- 
ulation responsible for this reionization remains an open 
question. The evolution in the bulk contribution from 
LAEs has only recently been explored and shown to be a 
substantial contributor of ionizing flux from z = 3.1 — 5.7 
(O08; M08). However, such samples are only able to di- 



rectly measure contributions from galaxies on the bright 
end of the luminosity function. The contribution from 
sub-L, LAEs are inferred by extrapolation of the best-fit 
luminosity function to faint Lya line luminosities. Since 
many of our candidates lie far below the typical esti- 
mates for at these epochs, we can better characterize 
the contribution from such populations. It may be that 
galaxies with typical Lya luminosities comparable to our 
sample (< O.li* — L^) represent the dominant contribu- 
tion to the cosmic SFR amongst LAEs, a trend observed 
in z=2-5 LBG populations (Sawicki & Thompson 2006). 

Using the entire survey volume and all our LAE can- 
didates, we find an SFRD of 5.2;j ° x lO'^ (4.5;5j;^ x 

10~^ excluding Q=l candidates) Mq yr~^ Mpc"'^. These 
values should be viewed as lower limits to the SFRD, as 
we make no corrections for added slit loss (due to the un- 
known position of the LAE candidate) or correction for 
extinction of Lyman continuum photons. Though dust 
corrections are important to determine the SFRD, this 
correction is perhaps not problematic if we are concerned 
only with the number of photons available to ionize the 
universe at these epochs as any Lya photon that is ab- 
sorbed by dust will make little or no contribution to the 
re-ionization of the universe. The errors should also be 
viewed as lower limits, as we do not incorporate formal 
errors due to cosmic variance, though we make some ef- 
fort to quantify its effects (see Section 5.3). Despite these 
complications, a comparison the observed SFRDs of the 
C11604 LAE candidates gives values consistent with or 
exceeding the contributions of super-L, LAE galaxies 
found a,t z = 5.7 (M08, O08), but significantly less than 
the contributions from z = 3 — 6 LBGs (~ 1 — 2 x 10^^ 
M© yr~^ Mpc~^: Giavahsco et al. 2004; Bouwens et 
al. 2004a, 2006; Sawicki & Thompson 2006, Iwata et al. 
2007). While some LBG surveys correct for contributions 
from galaxies dimmer than the completeness limit of the 
survey by integrating the observed luminosity function, 
we have made no such correction. Even though this cor- 
rection could, in principle, be large enough to push our 
observed SFRD to levels competitive with LBG surveys, 
the faint-end of the LAE luminosity function is less con- 
strained than the low luminosity end of the LBG lumi- 
nosity function making extrapolation uncertain. If we 
assume that LAEs behave like LBGs at low luminosities, 
exhibiting relatively shallow faint-end slopes, the con- 
tributions from galaxies dimmer than the completeness 
limit of our survey (~0. 1-0.2 L^) would contribute only 
10%-15% to the total SFRD at the redshifts of interest. 

To quantify the evolution across the redshift range of 
our LAE candidates, we split our data into two redshift 
bins dictated by two OH transmission windows: (1) from 
z = 4.1 (the onset of our spectral sensitivity) to z = 
4.95 (the onset of significant contamination from bright 
airglow lines, see Figure [1]), and (2) between z — 5.6 
and z = 5.8, an atmospheric transmission window used 
by many narrowband imaging surveys. These choices of 
bins exclude only one LAE candidate, 16XR1.97 at a 
redshift of z = 5.02. 

The first redshift bin contains 14 LAE candidates, 11 
of which are high quality. The volume of the survey in 
this wavelength range, calculated in the same manner as 
in Section 2.2, is 4.99xl03 Mpc^ The SFRD for the 
lower redshift sample is SFRDz^4.55 = 12. it^;^ x 10~^ 
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{ll.ltlj X 10"^ excluding Q=l candidates) M© yr^^ 
Mpc~'^, a density rivaling the contribution of LBGs at 
this epoch. While this bin contains a large fraction of 
our LAE sample, any conclusions must be tentative, as 
cosmic variance can dramatically change the observed 
value (see Section 5.3). 

The higher redshift bin contains two candidates (both 
high quality) within a survey volume of 1.44x10^ Mpc"^. 
The SFRD for the higher redshift sample is SFRDz^s.y 
= 4.4t}'i X 10~^ M0 yr~^ Mpc~^, more consistent with 
the average SFRD of the survey and consistent within 
the errors of extrapolated SFRDs found by other sur- 
veys at similar redshifts (Rhoads et al. 2003, M08, O08). 
Since the high redshift bin contains only two LAE can- 
didates, the measured value of the SFRD in this bin is 
highly susceptible to cosmic variance effects. While any 
conclusions that involve the high redshift bin are very 
uncertain, the drop in SFRD density at high redshift is 
a statistically significant effect and could possibly rep- 
resent real evolution in the properties of LAEs as the 
observational epoch nears the epoch of re-ionization. If 
we instead choose to exclude the high redshift candidates 
from our sample and only use the lower redshift bin at 
z ~ 4.55, the observed SFRD is still significantly larger 
than those of higher redshift samples of LAEs (Rhoads 
et al. 2003; Ajiki et al. 2003; Shimasaku et al. 2006, 
hereafter S06; O08; M08). However, the derived SFRD 
at z ~ 4.55 is also quite a bit higher than some samples 
at comparable redshifts (Dawson et al. 2007, S08), sug- 
gesting that the observed change in SFRD from z ~ 4.55 
to z ^ 5.7 probably arises through some combination of 
cosmic variance effects and real evolution in the LAE 
population. While marginally inconsistent with other 
measurements of the evolution of the LAE SFRD from 
z - 4.55 to z r^6 (03; 08; Rhoads et al. 2003; Dawson 
et al. 2007) , these results are consistent with the general 
evolution of the star-forming properties of LEG popula- 
tions (Bouwens et al. 2004a; Sawicki & Thompson 2006; 
Yoshida et al. 2006; Iwata et al. 2007) and the overall 
picture of decreasing contribution to the cosmic SFRD 
from LAEs with increasing lookback time (Taniguchi et 
al. 2005). 

4.4. Velocity Profiles 

The observed line profile of the unsmoothcd composite 
spectrum of the 13 high quality Lya emitters was fit us- 
ing a five parameter single Gaussian model similar to the 
model used in S08. In this model we assumed that the 
ISM absorbed all Lya photons blueward of the centroid 
of an unattenuated Gaussian emission line, allowing us 
to produce the characteristic shape of the Lya line. The 
mean wavelength of the original Lya emission was al- 
lowed to freely vary, as our uncertainty in the redshift 
is coupled to our inability to quantify the extent of the 
absorption blueward of the Lya line. Additionally the ef- 
fective FWHM of our spectral setup, in principle a known 
quantity, was allowed to vary due to our ignorance of the 
placement of the LAE on the slit and the magnitude of 
the change in FWHM resolution as the LAE moves out 
of the slit. The background, dispersion, and amplitude 
of the Gaussian were also allowed to vary. This is, of 
course, a very simple model of the Lya emission. In real 
galaxies there are typically multiple emission components 



offset in velocity space. In the case of LAEs, there can 
also be a significant excess of flux in the far red end of 
the line profile due to backscattering of Lya photons by 
surrounding H II regions (S08; Westra et al. 2005). StiU, 
this model allows us to gain some insight into the av- 
erage properties of the main velocity component of our 
LAE candidates. 

Figure [Tl] shows the best-fit model line profile overplot- 
ted on the co-added spectrum of the high quality LAEs. 
This simple model does reasonably well reproducing the 
observed line profile. It is interesting to note that if the 
model represents, even roughly, the intrinsic, unatten- 
uated Lya line, the truncation of the Lya line by the 
IGM results in an attenuated line which is offset from the 
original line profile by nearly 100 km s~^. The best-fit 
intrinsic velocity dispersion of 136 km s~^ is marginally 
consistent with the findings of S08 and LAEs detected 
in some narrowband imaging surveys (H04) and is at the 
extreme low end of the mass function of other surveys 
(M08; Dawson et al. 2004). 

There are two main discrepancies between the data and 
the simple truncated Gaussian model. The first is the 
failing of the model to reproduce flux just blueward of the 
centroid of the Lya line at ~ 1215 A and again at 1214 
A. Such excesses could arise from a non-trivial amount of 
Lya photons escaping attenuation from pockets of neu- 
tral hydrogen. Indeed, even at the highest redshifts of 
our Q=3 LAE candidates (z ~ 5.6), Equation [7] predicts 
an escape fraction (fesc) of ~ 5%, increasing to more 
than 40% at the lowest redshifts. The model also fails to 
produce enough flux at the extreme red end of the line, 
showing a moderately significant decrement in flux at ^ 
1217.5 A as compared to the data. This unaccounted 
flux could be explained by backscattering of Lya pho- 
tons from galactic outflows as a result of star formation 
processes (Dawson et al. 2002; Mas-Hesse et al. 2003; 
Ahn 2004; Westra et al. 2005, 2006; Hansen & Oh 2006; 
K06). The offset between the observed flux excess and 
the centroid of the Lya emission is ~440 km s~^, con- 
sistent with this interpretation and with measurements 
from other surveys (Dawson et al. 2002: 320 km s~^; 
Westra et al. 2005: 405 km s'^ S08: 420 km s'^). It is 
interesting to note that these signatures appear in both 
the luminosity and uniform weighted stacked spectra, 
suggesting that such outflow processes are pervasive in 
low-mass high-redshift star-forming galaxies. However, 
both excesses are near the level of the noise in the co- 
added spectrum. While it is plausible to attribute these 
excesses to such astrophysical processes, more data are 
necessary to make any definitive conclusions. We, there- 
fore, defer more complicated modeling of the composite 
emission line profile until all ORELSE fields are included. 

5. LYa EMITTER NUMBER COUNTS AND LUMINOSITY 
FUNCTIONS 

5.1. Number Counts and Cumulative Number Density 

The depth of our data allows us to detect galaxies down 
to a limiting luminosity of 8.8 x 10^^ ergs s~^ over the 
entire redshift range of our survey and down to 3.7 x 10^^ 
ergs s~^ at z = 4.4. This is almost a factor of 2 deeper 
than the recent spectroscopic survey of S08, previously 
the deepest survey for LAEs to date, and nearly a fac- 
tor of 10 deeper than the completeness limits of recent 
narrowband imaging surveys (O08; Dawson et al. 2007; 



16 



velocity (km/s) 

500 



Original Gaussian 



: Truncated Smoothed 
Gaussian 



DEIMOS FWHM = 2.71 A 
= 1215.50 A 



ij a, = 135.75 (km/s) 




J Average 1 a Error 



1214 1216 
(A) 



io-« 



M0-» 



V 

E 



io-< 



Ouchi03, ; 
SawickiOB, <z>=4.85 ■ 
ShimasakuOe, z=5.7 - 
DawsonO?, z=4.5 
MartinOB, z=5.7 
KashikawaOe. z=6.5 ■ 



io-» 




10« 1043 

L(Lya) [ergs/s] 



Fig. 11. — A simple single Gaussian model fit to the unsmoothed 
composite (luminosity weighted) spectrum of high-quality (Q=2,3) 
LAE candidates. All flux blueward of the peak of the Gaussian 
has been removed in order to approximate the effects of attenu- 
ation from HI regions. The resultant proflle is smoothed with a 
second Gaussian simulating instrumental broadening. The peak 
wavelength and width of the original Gaussian as well as the in- 
strumental broadening are all free parameters in the model. The 
model proflle fits the data extremely well, with a few notable ex- 
ceptions. 

Murayama et al. 2007). Since our survey probes deeper 
in the luminosity function than previous surveys, any re- 
sults that involve raw number counts of LAE candidates 
must be corrected for differing survey flux limits if proper 
comparisons are to be made. 

One way to disentangle number counts from the ef- 
fects of varying flux limits is to consider the cumulative 
number density of LAEs. Since the overall shape of the 
function should be identical, in the absence of cosmic 
variance and any instrumental effects, the galaxy popu- 
lations from surveys of differing flux limits can be cast 
in a single functional form. Figure [12] shows both the 
cumulative number counts of all LAE candidates and of 
the highest quality (Q=3) LAE candidates in the C11604 
field as compared to other surveys. Although the bulk 
of the candidates in the C11604 field reside at z < 5, the 
number counts lie above the measurements of surveys at 
similar redshifts. These number counts are also signifi- 
cantly higher than the lower limits of SOS, a survey with a 
nearly identical instrumental setup to our own. Extrapo- 
lating down to the limiting luminosity of our survey, the 
C11604 LAE candidates seem to be instead broadly con- 
sistent with the number counts of LAE surveys at much 
higher redshifts. 

There are many possible explanations for this dis- 
crepancy. While more exotic possibilities are discussed 
in Sections 5.4 and 5.5, there are simpler explanations 
which can be explored immediately. One possibility that 
would make our results slightly more consistent with 
other surveys at similar redshifts is the exclusion of the 
object 16XR1.72 (the LAE candidate at LLj^a~l-7x 10''^ 
in Figure fT2| . This is an extremely bright (~2-3 L*) ob- 
ject that is unlikely to have been detected in a survey 
with our limited volume. Its presence is a strong indi- 
cation that cosmic variance of such galaxies may have 



Fig. 12. — Cumulative number density of LAEs detected in 
the C11604 field plotted as a function of hya line luminosity mini- 
mally corrected for slit losses {u> = 0.8; see Section 2.3). The right 
axis shows the cumulative number of LAE candidates detected in 
our survey. The observed number counts are more consistent with 
LAE populations at higher redshift (z > 5.7) than populations at 
more moderate redshifts (z < 4.9). Even excluding all but our 
highest quaUty (Q=3) candidates the number density is consistent 
with only the highest measurements of other surveys at similar red- 
shifts. Error bars are derived from a combination of Monte Carlo 
simulations that incorporate the uncertainties in the luminosities 
(effectively which bin a given LAE candidate will fall in) and Pois- 
son statistics, assuming (improperly, see Section 5.3) that LAEs 
have a spatial distribution reasonably similar to that of a Gaussian 
random field. Introducing formal measurements of uncertainties 
due to cosmic variance (not yet possible for this sample) will cause 
these errors to increase. Note that the SOS number counts plotted 
in this figure and Figure [141 differ from those plotted in SOS as we 
include LAEs detected at all redshifts in their survey (with the 
exception of two marginal candidates), whereas they only include 
LAEs detected between z = 4.2 and 4.9 

biased our results high for bright LAEs and its exclusion 
would serve to eliminate this effect at the brightest end 
of the cumulative number density distribution. Another 
possibility is the large uncertainties in the fluxes (lumi- 
nosities) of the LAE candidates due to slit loss and flux 
calibration. As discussed in Section 2.3, the absolute flux 
calibration of the data was accurate to only 60%, further 
compounded by a possible systematic offset which un- 
derestimated the true value of the flux. Ignoring for the 
moment the systematic offset, the uncertainty in the flux 
calibration, combined with Poisson errors, could cause 
any given LAE candidate to shift nearly a factor of 2 in 
brightness (Icr) in either direction. Although significant, 
this shift would not change our conclusions, as our data 
would still be broadly consistent with the higher redshift 
surveys and broadly inconsistent with the moderate red- 
shift surveys. Any systematic offset due to slit attenua- 
tion (see Section 2.3 for a full analysis) would cause an 
underestimate in the line fiux and shift the cumulative 
number density curve to the right, only reinforcing our 
conclusions. If flux errors are the only potential contam- 
inant in our results, our data suggest minimal evolution 
in the Lya number density from z — 4.4 to z = 6.5. How- 
ever, as we explore in the following two sections, there 
are other explanations which allow for evolution in the 
hya luminosity function. 
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Fig. 13. — The LAE candidates plotted against the foreground 
of the C11604 supercluster, with the two possible LAE structures 
at z ~ 4.4 and z ~ 4.8 shown as separate symbols. A redshift 
histogram of all LAE candidates (more finely gridded than the 
one in Figure (7)1 is shown to the right of the supercluster. For 
illustration, the name of each cluster as well as dashed circles that 
represent radii of 0.5 h~^ Mpc are overplotted. The dot-dashed and 
dotted lines in the histogram show the redshift distributions of the 
possible LAE structures at 2 ~ 4.4 and z ~ 4.8, respectively. The 
solid line shows the redshift distribution of the remaining "field" 
LAE candidates. 



5.2. Clustering of Lya Emitters at z ^ 4.4 and z ^ 4.83 

Since our survey encompasses an extremely small vol- 
ume compared to other contemporary surveys, we are 
susceptible to sampling galaxies whose distributions are 
unrepresentative of the distributions observed in larger 
LAE surveys. Effects such as strong clustering of LAEs 
(O03) which lead to high levels of cosmic variance can 
play a large role in surveys with limited breadth. Simi- 
larly, variance in the observed number densities of LAE 
populations could be caused by inhomogeneities of inter- 
vening HI regions, as areas with a sparser density could 
manifest in increased detections of LAEs. While we will 
defer any complex treatment of cosmic variance until we 
can include all of the ORELSE fields, we attempt to 
quantify its effects in our survey in the following sec- 
tion. However, preliminary results from LAE searches in 
other ORELSE fields indicate similar detection frequen- 
cies, suggesting that cosmic variance may not be the sole 
cause of the large LAE number densities observed in the 
C11604 field relative to other z ~5 surveys. 

We do see significant evidence for clustering in our 
data, so it is possible that this clustering, combined with 
strong (Poissonian) sample variance, could be enough to 
explain the observed discrepancy in LAE number counts 
in the C11604 field relative to other surveys at the same 
redshift. As shown in Figure [71 there is a clear redshift 
peak at z ~ 4.8, as well as a less pronounced peak at 
z ~ 4.4. This is recast in Figure [13] where each LAE is 
plotted against the foreground of the supercluster, with 
the two redshift peaks being differentiated from the gen- 
eral LAE population. 

While the galaxies are not strongly clustered spatially, 
the observed distribution is broadly consistent with the 
spatial distributions and number densities of other struc- 



tures found at high redshift (e.g., Shimasaku et al. 2003) 
and represents a significant overdensity when contrasted 
with normal "field" populations of LAEs. Using the cov- 
erage provided by the slits to calculate the survey vol- 
ume at these redshifts, we recover a cumulative number 
density of 2.6 x 10"^ Mpc"^ and 1.4 x 10"^ Mpc~^ for 
the z ~ 4.4 and z ~ 4.83 structures, respectively, which 
are inconsistent at the >99.99% C.L. with field counts of 
LAEs (see Figure (Til) . 

Since we have truncated the bounds of the data know- 
ing the redshift range of the two structures involved, it is 
possible that we are overestimating the number density 
for these structures. However, in the absence of clus- 
tering we would expect our data to closely resemble the 
DEIMOS response function in areas absent of bright air- 
glow lines. Since we are sensitive to a level far below 
typical estimates of L* across the entire redshift range of 
the survey, there should be little dependence of the num- 
ber counts in our data on epoch. We are most sensitive to 
the redshift range 4.5-4.9 where the DEIMOS through- 
put is maximized, the airglow lines are minimal, and the 
full survey area is exposed to these wavelengths. It is 
surprising, therefore, that nearly all of our candidates in 
this redshift range lie in a narrow peak in redshift space, 
and that another peak exists outside the range of this 
area of high spectroscopic sensitivity. 

If we exclude the LAEs in these two redshift peaks from 
our general LAE sample on the grounds that they belong 
to rare structures, we are left with a total of seven (six 
high quality) LAE candidates in our sample that rep- 
resent a typical sampling of the LAE field populations. 
This reduction in the number of candidates drops our 
number counts to levels consistent with, but still higher 
than, other surveys at similar redshifts within the bounds 
of reasonable (Poissonian) sample variance (see Figure 

m- 

5.3. Cosmic Variance 

Since the LAE candidates detected in the C11604 field 
represent one of the first detections of a reasonably large 
sample of faint LAEs, little is known about the cluster- 
ing behavior of such galaxies. While the similar depth 
of the spectroscopic data in other ORELSE fields will 
allow for the study of statistical properties of faint LAE 
populations, the data presented in this paper is limited 
to one field. Thus, it is possible, given the clustering 
observed of brighter LAEs (O03; Shimasaku et al. 2003; 
S06), that our measurements of cumulative number den- 
sities (Section 5.1) or luminosity functions (Section 5.4) 
are biased by uneven sampling of the spatial distribution 
and the dumpiness of the population. As large statistical 
samples of faint LAEs do not exist, we must estimate by 
other means the likelihood that we would have detected 
the same number of LAE candidates if our survey had 
observed a different region of the cosmos. 

In order to estimate this likelihood and to measure 
the magnitude of cosmic variance on the C11604 LAE 
candidates, simulations were run on four different sam- 
ples of LAEs. The four fields were chosen because they 
contain a large number of LAEs, which were uniformly 
(or nearly uniformly) sampled over a large comoving vol- 
ume and spanned the redshift range of the C11604 LAE 
candidates. The samples were: 1) the LALA spectro- 
scopic sample in the Bootes field at z ^4.5 (D07), 2) 
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F ig. 14. — The cumulative number density of LAEs as in Fig- 
ure [TS] but with the two possible LAE structures at 2 ^ 4.4 and 
z ~ 4.8 differentiated from the general population. The volume 
for each potential structure is bounded by the redshift range of the 
LAE candidates in the sample and calculated from the slit area 
exposed at those redshifts (wavelengths). The number densities of 
these potential structures are several orders of magnitude above 
field populations at all redshifts. If these structures are real, the 
number density of our remaining LAE candidates begin to be more 
consistent with moderate redshift (2 < 4.9) LA E po pulations. Er- 
ror bars are derived in the same way as Figure 1121 Error bars in 
the last bins of the two LAE structure curves have been removed 
for clarity. 

the Subaru Deep Field (SDF) at z -4.9 (03; hereafter 
SDF+LSS), 3) the SDF at z ~4.9, excluding the vol- 
ume containing the large scale structure as defined in 
Shimasaku et al. (2003; hereafter SDF-LSS), and 4) the 
SDF at z -5.7 (S06; hereafter SDFhighz). The LAE 
candidates in all fields were selected using narrowband 
imaging methods and, in some cases, followed up with 
spectroscopy. Although each survey sampled a large vol- 
ume (~lxl0^ Mpc'^), their coverage was concentrated 
at nearly one epoch. However, since each survey has 
moderately large coverage in the transverse directions, 
the variance observed in each sample likely represents a 
reasonable estimate for the cosmic variance of brighter 
LAEs at that epoch. For simplicity, the transmission 
of the narrowband filter was assumed to be a top-hat 
response, with the hat size equal to the FWHM of the 
true filter response curve centered around the effective 
wavelength. 

For each realization of the simulation in each field, an 
area was "observed" that would yield the survey volume 
of the C11604 spectral data (i.e. 1.365x10'' Mpc^) given 
the filter setup. For the SDF samples, where LAE can- 
didates were selected using only one narrowband filter, 
this area corresponded to roughly 65 arcmin^. In the 
LALA field the area observed in each simulation was 
significantly less (— 22 arcmin^) due to the continuous 
coverage of their five narrowband filters, which span the 
line of sight direction from z — 4.37 to z = 4.57. The 
observation in each realization consisted of counting the 
number of LAEs detected in a continuous square area, 
whose central R.A. and decl. were determined by ran- 
domly drawing from a uniform distribution bounded by 
the spatial coverage of each survey. The results of these 



simulated observations are shown in Figure [TSl 

The results of each simulation are considerably dif- 
ferent, suggesting that, even in surveys for LAEs that 
probe large comoving volumes, cosmic variance can play 
a large role or, alternatively, that the clustering statistics 
of LAEs evolve rapidly between z —4.5 and z —5.7. The 
variations may also arise from the difference in the pa- 
rameters of each survey (e.g., limiting magnitude, com- 
pleteness, purity). The results using the LALA sample in 
the Bootes field likely represents a lower limit to the cos- 
mic variance and simulated LAE number counts because 
only 60% of LAE candidates selected by narrowband 
imaging were targeted by spectroscopy. Conversely, both 
the SDF+LSS and SDF-LSS samples have estimated pu- 
rities of 60%-70% (O03). Thus, they represent an up- 
per limit to the cosmic variance. The SDFhighz sample, 
which consists of a mix of spectroscopically confirmed 
LAEs and objects selected solely through narrowband 
imaging, falls somewhere in between the z — 4.9 SDF 
and the LALA sample in terms of completeness and pu- 
rity. 

In each panel of Figure [151 the two dashed lines mark 
the number of LAEs that we expect to observe in the 
C11604 field (i.e., 3) based on the extrapolating the num- 
ber counts of surveys of LAEs at similar epochs (D07; 
SOS) to the limiting luminosity of our data and the actual 
number of LAE candidates (i.e., 17) that we detect of all 
qualities. For both the LALA and SDFhighz fields, these 
simulations rule out cosmic variance as the sole cause of 
the observed excess of LAEs detected in the C11604 field 
at > 99.99% C.L. However, since the bulk of our galaxies 
lie between z—AA and ^=4.9, and since the LALA data 
are sparsely sampled, we focus on the results of the two 
SDF fields at z —4.9. In both cases, the observed num- 
ber of LAEs detected in the C11604 are allowable within 
the bounds of the simulated cosmic variance. The like- 
hhood of recovering at least 17 LAEs in the SDF-fLSS 
simulations is 26% compared to only 6% in the SDF-LSS 
simulations. This result strongly supports the conclusion 
reached in §5.2: we may be observing at least one large 
scale structure of LAEs in the C11604 field. 

However, each LAE sample from which these simula- 
tions were drawn is, on average, significantly brighter 
than the C11604 LAE candidate population. The com- 
pleteness limits of each survey are roughly 3-8 times 
brighter than that of the C11604 survey. Thus, many 
of the C11604 LAE candidates may not be detected in 
these surveys. In order to place a lower bound on the 
number of galaxies that would have been detected (as- 
suming all of our LAE candidates were at the survey 
redshift), we cut the C11604 LAE candidate population 
at Lya luminosities at or above the luminosity corre- 
sponding to the completeness limit of each survey. This 
number, plotted as a dotted line in each panel of Figure 
[TSl is a lower limit since the luminosities calculated for 
the C11604 LAE candidates are lower limits. Since we are 
cutting at the completeness limit and not the limiting lu- 
minosity, the number of LAEs that would be detected by 
each survey had they observed the C11604 field and cov- 
ered a volume equivalent to the C11604 survey volume lies 
somewhere between the dotted and the rightmost dashed 
line in each panel. Including this cut, we find that the 
hypothesis that our field contains a large-scale structure 
is still favored, though less strongly, as 44% of the sim- 
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Fig. 15. — Four different simulations of the effects of cosmic vari- 
ance on the observations in the C11604 field. The four simulations 
are designed to cover a large range of redshifts and LAE samples. 
Each panel represents simulated observations of LAEs from differ- 
ent narrowband imaging surveys (see Section 5.3 for details on the 
surveys) . A histogram of the number of LAEs recovered in each re- 
alization is plotted in each panel. The two dashed lines correspond 
to the number of LAEs expected in the C11604 data sampling an 
"average" field at z ~5 (Nexp) and the number of actual detections 
(all qualities) of LAE candidates in the C11604 data (Nobs)- The 
dotted line in each panel represent the number of LAEs that we 
would have detected in the C11604 field if we instead adopt the 
completeness limit of each survey. These results suggest that cos- 
mic variance is a major contributor to the excess of LAEs described 
in §5.1. 

ulations in the SDF-I-LSS field recover 12 or more LAEs 
(the number of C11604 LAE candidates above the O03 
luminosity limit) compared to only 25% in the SDF-LSS 
field. 

While it is unclear what adding in fainter LAEs 
to these simulations would do to the measurement of 
the cosmic variance, it has been observed that low- 
luminosity (~0.3-0.5L,) LAEs are less strongly clustered 
than brighter (> 0.5L*) LAEs (O03). While this trend 
may not extrapolate down to the limiting luminosities 
of this survey or may be an effect unique to the SDF 
at z ~4.9, the consequences of adding in the population 
may be limited and may even serve to dilute cosmic vari- 
ance. While no definitive conclusions can be reached, the 
main result of these simulations is that cosmic variance 
may be solely responsible for the observed excess of LAE 
number density detected in the C11604 field. Determin- 
ing how likely that is, however, is beyond the ability of 
these simulations. 

5.4. Lya Emitter Luminosity Function at (z) — 4.85 

Because a large majority of our LAE candidates re- 
main undetected in our imaging, these objects are equally 
likely to fall anywhere within the bounds of our slits 
or perhaps, depending on their brightness, outside the 
bounds of our slits. Without any way of recovering the 
true position of these objects, the slit loss correction must 
be approached statistically. 

The statistical correction is made using the simula- 
tion code designed specifically for this purpose by Martin 
and Sawicki (see MS04 for a more detailed explanation). 
Briefly, for each realization the underlying galaxy pop- 



ulation is simulated by setting the parameters L*, 
and a, fully describing a unique instance of the under- 
lying distribution characterized by the Schechter (1976) 
function: 



exp 
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dL 
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Data are then simulated for a grid of Schechter param- 
eters for redshift slices of 6z — 0.1 running from central 
redshifts of z = 4.15 to z = 6.35 and multiplied by the 
area exposed to each redshift interval in our survey. The 
area exposed on the sky was also allowed to vary (in a 
known way) as a function of source flux and the seeing. 
Since brighter objects can fall further from the center of 
the slit (widthwise) and still be detected by our survey, 
the area of our survey at all redshift intervals increased 
with increasing LAE flux. The range of simulated LAE 
luminosities for each redshift bin was left unbounded on 
the bright end and truncated on the faint end by an LAE 
luminosity that would result in a flux of 1.9 x 10^^^ ergs 
s~^ cm^^ (ranging from L(LyQ;) = 3.2 x 10'*^ ergs s^^ at 
z = 4.15 to L(Lya) = 8.9 x lO^i ergs s"! at z = 6.45). 
For every set of Schechter parameters, each simulated 
galaxy is "observed" by calculating the slit attenuation 
based on simulated slit losses for an unresolved source 
galaxy at a regularly sampled grid of positions with re- 
spect to the slit in 0.9" seeing. The total number of 
LAE galaxies of all fluxes (luminosities) for that set of 
Schechter parameters is then recorded. Though we only 
include galaxies in these simulations with fluxes greater 
than or equal to the completeness limit calculated in Sec- 
tion 2.5, this choice results in conservative values of the 
Schechter parameters. Specifically, we underestimate the 
"true" Schechter parameters since the actual complete- 
ness limit of our survey is brighter than the limit calcu- 
lated in Section 2.5 (due to the unknown amount of slit 
losses). Therefore, it is likely that the faintest galaxies 
in this simulation were not detected in our survey. Thus, 
this choice results in more simulated galaxies than if we 
had attempted to make a correction for slit losses, re- 
quiring us to observe more galaxies in the C11604 field to 
recover the same set of Schechter parameters. 

From this simulation we are not able to recover a 
unique set of Schechter parameters due to our ignorance 
of how many genuine LAE galaxies are in our data. Fur- 
thermore, since the simulations allow each LAE candi- 
date to have a range of Lya luminosities, we are not 
definitively setting the number of galaxies detected at 
any given luminosity. This constraint would be essen- 
tial if we were calculating specific values of and L*. 
Instead, we can only limit the range of Schechter param- 
eters by bounding the number of simulated galaxies for 
a given ^^L^ by the total number of genuine LAE galax- 
ies in our data. It seems reasonable that the number of 
galaxies observed in the simulation be equal to at least 
13, the number of high quality LAE candidates in our 
sample. However, this number does not represent a hard 
lower bound, as any real clustering in the data would 
not be accounted for in these simulations. Instead we set 
the lower bound as seven (six high quality) galaxies, the 
number of LAE candidates that exist outside the two 
possible structures and constitute a lower limit on our 
LAE field population. 
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Fig. 16. — The range of possible Schechter parameters, 
and L«, for simulated populations of LAEs approximating those 
observed in the C11604 field for two different faint-end slopes 
(a = —1.2, top panel; a = —1.6, bottom panel). The data are 
simulated assuming a completeness limit of 1.9 X 10~^* ergs s~^ 
and is corrected for statistical flux losses. Each contour through the 
^tLt phase space represents the expected number of LAEs, given 
our instrumental setup and observing conditions, that should be 
detected in our data. The shaded contour shows the phase space 
allowed by our LAE candidates. For comparison the values of 
and L, are shown for different surveys that used similar values 
of the faint-end slope. (MR04: Malhotra & Rhoads 2004; K06: 
Kashikawa et al. 2006; S06: Shimasaku et al. 2006; D07: Dawson 
et al. 2007; G07: Gronwall et al. 2007; O08: Ouchi 2008) 

The upper bound for this simulation is more ambigu- 
ous. The primary consideration is the nature of our com- 
pleteness limit, the flux at which galaxies were cut in the 
simulation. We estimate that at our flux limit we do not 
miss more than two-thirds of the actual number of LAE 
galaxies in our data. Taking the 17 LAE candidates of all 
qualities as the upper bound on actual detections, this 
sets an upper bound for our simulations at 51 LAEs. 

For each panel in Figure fT6l the shaded area represents 
the values of the Schechter parameters, the normalization 
and the characteristic luminosity (L*) allowed by 
our data for the range of possible LAE galaxies detected 
by our survey. This analysis is done for two different faint 
end slopes, a=-1.2 and a=-1.6. Even if we knew the true 
number of LAEs in our data, due to the nature of the 
Schechter formalism and the sparseness of our data, we 
are not able to definitively determine unique values of the 
Schechter parameters from our simulations. However, 



our data does allow us to constrain the product, <1>*_L*, 
by the range bounded by the shaded area in Figure [TH 
Defining the best-fit Schechter product as the average 
value of for which we observe exactly 13 LAEs in 

the simulation and the range of possible •I'^L* values as 
the average of those which recover seven simulated LAEs 
(on the low end) and 51 simulated LAEs (on the high 
end), we find a best-fit Schechter product of $*L* — 
2.2^11 X 10^^ ergs s~^ Mpc~^ for a faint-end slope of 
a = —1.6. 

As shown in Figure \W[ our observed range in and 
L^, is consistent with measurements made by D07, O03, 
and MR04 and slightly low when compared to measure- 
ments made by K06, S06, G07, and O08. Although there 
is considerable variation in the measured values of 
and even at similar redshifts, the comparisons of the 
Schechter parameters seem inconsistent with the com- 
parisons in Section 5.1, in which the cumulative number 
density of the C11604 LAE candidates were more simi- 
lar to the extrapolated number counts of higher redshift 
surveys (K06, S06) than surveys at lower redshift (e.g., 
D07). It is possible that our conservative estimates of 
both the number and luminosity of our candidate LAEs 
may be the source of this discrepancy rather than any 
real evolution in the luminosity function of Lya emitters 
from z ~ 4.85 to z ~ 6. Since the choice of a faint-end 
cutoff can severely affect simulated numbers of LAEs, es- 
pecially for steeper faint-end slopes, our choice of a sim- 
ulated flux limit of 1.9 x 10~^* ergs s""'^ cm~^ could have 
biased our results to lower values of <I>*L*. To estimate 
the magnitude of this effect we ran the simulation again 
on a small portion of the data {z — 4.1 to z = 4.9) with 
a brighter completeness limit of 3 x 10^^* ergs s~^ cm^^. 
The recovered Schechter product was on average higher 
by a factor of ~ 2, corresponding to a « 5 x 10"^^ 

ergs s~^ Mpc~^ for a faint-end slope of a = —1.6, essen- 
tially pushing the contours up and to the right in both 
panels in Figure [TH encompassing the <I>*L* products 
of the other surveys within our range of allowed values. 
Since the true completeness limit of the C11604 spectral 
data is somewhat uncertain (see Section 2) and since the 
results of these simulations are extremely sensitive to the 
choice of this limit, we are not able to distinguish be- 
tween the luminosity function properties of our sample 
and other samples of LAEs. Instead, we conclude that 
the Schechter parameters for the C11604 LAE population 
are broadly consistent with other measured values from 
z ~ 4 to z ~ 6. 

5.5. Weak Lensing Contributions to the Luminosity 
Function 

The C11604 supercluster is the most well-studied large- 
scale structure at high redshift. While no single cluster 
in the structure would be considered at the high end of 
the cluster mass function (Poggianti et al. 2006, 2008; 
Milvang-Jensen et al. 2008; Hamana et al. 2008), with 
the possible exceptions of cluster A [a^ ^ 703 km s^^) 
and cluster B (cr ~ 800 km s~^) (G08), the large number 
of moderately massive constituent clusters and the struc- 
ture's large spatial extent make it an efficient astrophys- 
ical lens. The nature of the source population also lends 
itself to a large lensing effect: the lensing efficiency for 
LAEs being a monotonically increasing function of red- 
shift for ZLya > ziens- The presence of this massive lens 
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along with the high lensing cfRciency for high redshift 
LAEs makes it necessary to properly account for lensing 
processes and determine whether such processes may ex- 
plain the observed excess in number density counts over 
comparable field studies (see Section 5.1). 

We consider the effects of strong and weak gravita- 
tional lensing separately, with strong gravitational lens- 
ing effects discussed in Section 5.5. The observational 
weak lensing effect considered here is a magnification of 
the source population due to the lensing-induccd increase 
in observed surface area while keeping surface brightness 
constant. Assuming the slit is sufficiently large (or equiv- 
alently the galaxy is sufficiently small) to encompass this 
increased surface area, this magnification increases the 
brightness of observed objects and the frequency of de- 
tection by: (1) by lensing objects into the slit which were 
not already within the detectable area of the slit and (2) 
by increasing the total flux of galaxies that were just 
below the detection limit of the survey. Weak gravita- 
tional lensing could then increase both the overall nor- 
malization, and the characteristic luminosity, L*, of 
the LAE luminosity function. These effects are opposite 
those of all other analyses and measurement techniques 
used in this study, which are intentionally designed and 
implemented to underestimate the line flux. Such effects, 
if significant compared to the other associated uncertain- 
ties in our measurements, could significantly alter our 
conclusions. 

In principle, the most accurate approach to quantify 
the weak lensing effect would be to correct each LAE for 
the lensing- induced magnification. However, much more 
spectroscopy in the field (or other similar data) would be 
necessary to accurately measure the effect of weak lens- 
ing on each LAE. While a formal weak lensing analysis 
has been done on a small subsection of the field around 
cluster A (Margoniner et al. 2005) and will be done again 
with newly obtained ACS data (Lagattuta et al. 2010), 
the current data require us to take a more general ap- 
proach. In order to simulate the effect of weak lensing, 
the eight clusters that comprise the C11604 supercluster 
were modeled by singular isothermal spheres (SIS) of the 
form: 

While this model is an oversimplification of the true 
cluster mass profile, an SIS was used in place of a 
Navarro-Frenk- White (NFW, Navarro et al. 1996) pro- 
file because of the closed analytic form of the convergence 
and magnification solution. More importantly, however, 
since some of the clusters in C11604 are poorly sampled 
by spectroscopy (clusters F, G, H, I, and J), we were un- 
able to constrain the characteristic radius, r,,, and the 
concentration parameter, c, needed to properly charac- 
terize an NFW profile. 

Each cluster profile was simulated using velocity dis- 
persions published in Gal et al. (2008) and central posi- 
tions determined from the velocity and spatial centroids 
of the constituent cluster members. A velocity disper- 
sion derived from cluster members within 1 Mpc was 
adopted, chosen over 0.5 Mpc or 1.5 Mpc be- 
cause it is the largest velocity dispersion that is relatively 
free from significant contamination from other clusters. 



Since p scales as cr^, this choice will allow us to measure 
the maximum possible (model dependent) lensing effect 
on the LAE population by the clusters. 

For each realization of the simulation, new velocity dis- 
persions were generated for each cluster by a Gaussian 
sampling of the published 1 Mpc velocity dispersion 
errors. Thus, the velocity dispersion of the ith cluster 
was given by: 

Crv,i = <^i,lMpc + f^i, (13) 

where rij is sampled from a Gaussian distribution with 
width equal to the velocity dispersion errors. This new 
velocity dispersion, along with the static central positions 
of each cluster, completely dictated the mass map for the 
field of each realization; any effects from substructure, 
other structures along the line of sight, or lensing due 
to individual cluster galaxies were completely ignored. 
The source plane was created by averaging the z = 3.7 
and z = 5.7 field LAE luminosity functions taken from 
the large sample from the SXDS (O08), as our survey 
marginalizes over any evolution in the LAE luminosity 
function. The resulting Schechter function is parameter- 
ized by = 5.55 X lO"'' Mpc-^\ = 8.5 x 10''^ with a 
faint end slope of a = —1.5. The source population was 
drawn from discrete luminosity bins of width 0.1 dex, 
which were evenly spaced between L ~ 10^^ Lq and L 
~ 10*^ Lq. These limits were chosen to span the entire 
observable range of luminosities in various surveys. The 
number of galaxies in each luminosity bin was given by: 

..£-r..(A)%^-)-, 

where F is the simulation volume of 3.85 x 10^ Mpc^. All 
galaxies in each bin were assigned the average luminos- 
ity in that bin, dictated by the bounds of the integration. 
Next, each galaxy was assigned a R.A. (a) and deck {5), 
generated randomly from the area bounded by (but not 
necessarily sampled by) the spectral coverage of the sur- 
vey. In addition, a redshift was assigned to each source 
galax;y, drawn in a uniform random manner in the red- 
shift range z = 4.4 to ^; = 6.5. In principle, we could 
have made the simulation more realistic by introducing 
evolution in the LAE population by using Schechter pa- 
rameters based on various surveys at various redshifts. 
However, since the difference between the lensed and un- 
lensed populations is relatively insensitive to our choice 
of Schechter parameters, we ignore this effect. 

Once the ^ 25000 generated LAEs were assigned 
unique redshifts, coordinates, and luminosities, each lu- 
minosity was converted into an observed flux using the 
luminosity distance of each LAE. With the source plane 
and lens plane completely constructed, the lensing for- 
malism could then be applied. For the SIS profile the 
shear and convergence are equivalent, given for the ith 
LAE as: 

7.. = = 0.9 (^) (15) 

where the )th index represents the induced shear or con- 
vergence from the jih cluster in C11604 and 9 is the an- 
gular separation of the ith LAE from the jth. cluster in 
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Fig. 17. — Observations of simulated populations of LAEs in the 
absence of any lensing effects (solid lines) and after being weakly 
lensed by the C11604 supercluster (dashed lines). The two sets of 
curves show the differential number as a function of Ly« line lumi- 
nosity prior to the flux cut (continuous and short dashed lines) and 
subsequent to the flux cut (long dashed and dot-dashed lines). For 
all line luminosities the weak lensing number counts are consistent 
with the unlensed population within Sit. The biggest difference 
between the number counts (~10%) occurs at the line luminosity 
corresponding to the flux limit of 1 X 10~^® ergs s~^ cm~^ and 
dropping to negligible differences (1-2%) at higher luminosities. 

arcseconds. The magnification of the iih LAE is then 
calculated by: 
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The flux of each LAE was then increased by its re- 
spective magnification. A flux limit of 1 x 10"'^* ergs s~^ 
cm~^, reasonably approximating the flux limit of our sur- 
vey, was imposed on both the original, unlensed source 
population and the newly generated weakly lensed pop- 
ulation. The lensed flux was re-converted to an appar- 
ent luminosity, with the resulting observed populations 
shown in Figure flTl 

The results conclusively demonstrate that the cluster- 
induced weak lensing effect is far too small to account for 
our increased number counts. The ensemble average in- 
crease in total detections from the unlensed to the lensed 
data is 6%, an effect which is consistent with unlensed 
number counts at the 2a level in most of the bins. The 
effect is small regardless of the luminosity of the lensed 
galaxy, with an average increase in number counts in 
each bin ranging between 2% for the brightest simulated 
galaxies and 8% near the characteristic luminosity, L*. 

Another way to quantify the magnitude of this effect 
is by the overall increase in the Schechter parameters 
and L*. While not intuitively obvious as a measure- 
ment in the overall increase (or decrease) in the number 
counts of LAEs at different luminosities, it will give us 
some insight into possible systematic errors (as a result 
of unquantifled weak lensing effects) in our flnal luminos- 
ity function parameters. A Schechter parameter model 
with a fixed faint end slope (a = —1.5) was fit to both 
the unlensed and lensed data using a minimization 
routine. As the two Schechter parameters, and L*, 



are degenerate with one another (increasing necessi- 
tates a decrease in the normalization in order to maintain 
constant numbers of galaxies), the quantity of interest in 
these fits is not the individual parameters but rather the 
product For the unlensed data the best-fit pa- 

rameters resulted in a $*L* = 4.0 x 10^^ ± 2.5 x 10^*, 
differing from the product of the original input Schechter 
parameters due to the way the simulation is coarsely 
binned. Errors were calculated from the covariance ma- 
trix and the difference between the input and measured 
Schechter parameters. The fit to the lensed data resulted 
in a $*L* = 4.3 x 10'^^ ±3.1 x 10^^, representing an over- 
all increase of 7.5%, but also consistent within the errors 
to the unlensed measurements. While a correction of 
this magnitude might be important for precision mea- 
surements, the other uncertainties in our data (e.g., fiux 
calibration or flux losses due to the slit) far outweigh any 
induced weak lensing signal. 

5.6. Strong Lensing Contributions to the Luminosity 
Function 

When a galaxy is strongly lensed, either by the clus- 
ter potential or by a massive foreground galaxy, multiple 
images of the background galaxy are created on the sky. 
Depending on the relative positions of the lensing po- 
tential and the source (background) galaxy as observed 
projected on the sky, the resulting images of the original 
galaxy can be either fainter or brighter than the original 
galaxy, thus changing the observed luminosity function. 
This effect can also serve to push galaxies that would 
otherwise be too faint to detect above the flux detection 
limit, an effect which has been exploited by several sur- 
veys attempting to detect galaxies at very high redshift 
(z = 6-10, Santos et al. 2004; Egami et al. 2005; Stark 
et al. 2007; Richard et al. 2009). It is also possible, 
though very unlikely, that the slit geometry is perfectly 
oriented to observe multiple images of the source galaxy, 
effectively increasing the frequency of LAE detections. 
While these effects may introduce severe bias when they 
occur, strong gravitational lensing has a comparatively 
small cross section, allowing us to explore the possible 
effect on the LAE candidate population in a much more 
direct way than the exploration of similar effects caused 
by weak gravitational lensing. 

For an SIS lens multiple images form only for galax- 
ies that lie within the Einstein radius {Oe)- All of the 
LAE candidates lie very safely outside of any reason- 
able estimate of the cluster strong lensing regime (see 
Figure [T5|) , meaning any cluster-induced lensing effects 
would be accounted for in the previous section's analysis. 
The galaxy that comes closest to being strongly lensed is 
16XR1.26, which falls greater than 39e outside of clus- 
ter B. While the true mass profile may differ from SIS 
and the cluster mass may be underestimated due to our 
choice of velocity dispersions, it is unlikely that either of 
these effects would be strong enough to increase 9e to 
encompass 16XR1.26. 

Another issue that should be considered when dis- 
cussing strong lensing of the LAE candidate population 
is any strong lensing due to individual foreground galax- 
ies. This search has, by definition, a foreground galaxy 
companion, typically quite massive, targeted by the spec- 
troscopy. Since this galaxy lies within the length of the 
small DEIMOS slit, it is possible that this effect could be 
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Fig. 18. — LAEs with cluster members plotted with cluster 
Einstein ring radii {9e), which characterizes the onset of the strong 
lensing regime. Einstein radii were calculated using an SIS profile, 
with velocity dispersions derived from all members within 1 
Mpc of the cluster center. All LAEs fall clearly outside the bounds 
of the cluster-galaxy strong lensing regime. The galaxy closest to 
any cluster Einstein radius is shown in the bottom panel, lying 
many Einstein ring radii outside the center of cluster B. 

significant. In some cases there is an additional serendip- 
itous detection of a foreground galaxy on the slit which 
further complicates the matter. However, there are two 
reasons why we can ignore this effect. The Einstein ra- 
dius for an galaxy is of the order of 1". While this 
value will change based on the mass of the spectroscopic 
target (or other foreground serendip) and the relative 
redshifts of the foreground galaxy and the LAE, it is 
a reasonable estimate of where the strong lensing effect 
might be significant. All of our highest quality (Q=3) 
LAE candidates lie outside the bounds of this cutoff; 
any galaxy that is considered to be highly likely to be 
a LAE but is within a 1" radius of either the target or 
a foreground scrcndipitously detected galaxy is demoted 
to a lower confidence class (Q=2). Even considering only 
Q=3 LAE candidate galaxies, the detection frequency in 
our survey is significantly higher than most other surveys 
at the redshift of our sample, an effect which cannot pos- 
sibly be attributed to lensing. Furthermore, even if there 
are galaxy-galaxy lensing effects for which this analysis 
has failed to account, our sample set is selected nearly 
identically to the LAE population detected in SOS, with 
the possible exception that the spectroscopic targets in 
the C11604 field may be slightly more massive than those 
targeted by the DEEP2 survey. As they do not see sim- 
ilar excesses in their data it is likely that our observed 
excess comes from some combination of real, inherent 
properties of the observed LAE population, such as those 
discussed in the previous sections, and cosmic variance 
effects and cannot be attributed solely to lensing effects. 

6. SUMMARY 

In this paper we have described a search for LAE galax- 
ies in the 3.214 arcmin^ ORELSE spectroscopic database 
in the C11604 supercluster field. In total, 17 high redshift 
candidate galaxies were found in a volume of 1.365 x 10"' 
Mpc'^, with 13 galaxies meeting our high quality crite- 



ria. The redshifts of our LAE candidates ranged from 
2 = 4.39 to 2; = 5.67. Many of our candidate galaxies 
('--^90%) are dim compared to the typical characteristic 
luminosity at z ~ 5, with Lya line luminosities rang- 
ing from 5.9xl0''i erg s'^ (- O.IL,) to 1.7xl0''3 erg 
s~^ (~ 2L*). We have contrasted our LAE candidates 
with a population of known low redshift single-emission 
line interlopers and blended [Oil] emitters at intermedi- 
ate redshifts. Our 13 high quality candidates have prop- 
erties that differ significantly from the interloper popu- 
lation, giving us confidence in these objects as genuine 
LAE galaxies. The four lower quality objects do not 
distinguish themselves as well from the interloper pop- 
ulation, implying that these galaxies probably represent 
a mixture of LAEs and lower redshift objects. The in- 
creased frequency in LAE detections compared to other 
surveys demonstrates the effectiveness of LAE searches 
that probe deep into the luminosity function rather than 
covering large comoving volumes. Our main results are 
as follows: 

Lower limits on the Lya EW have been derived for 
all of our LAE candidate galaxies, finding a distribution 
peaking at EW(Lya)~20 A, similar to other low lumi- 
nosity galaxies at high redshift. We have also derived a 
lower limit to the SFRs of our LAE candidate, finding 
that they typically form stars at a rate of 2-5 Af© yr"^. 

From the entirety of our sample we determine an SFRD 
of 4.5tQ g X 10-3 Mq yr-i Mpc'^. This density is sim- 
ilar to or exceeding the contribution from super-L* LAE 
galaxies found at comparable redshifts, suggesting that 
sub-L, LAEs play an important role in keeping the uni- 
verse ionized at z 5. Grouping our LAE candidates 
into low redshift (4.1 < z < 4.95) and high redshift 
(5.6 < z < 5.8) bins, we find moderate evidence for neg- 
ative evolution in the SFRD. We measure an SFRD of 
2.6 1 n-3 at z 4.65 M© yr'^ Mpc'^ decreasing to 

M© yr^-' Mpc"'^ at z '--^ 5.7, though our 
highest redshift bin contains only two galaxies making 
any conclusions about the evolution of the LAE SFRD 
tentative. The derived SFRD of LAEs at z 4.55 is 
nearly equivalent to contributions of LBGs at similar 
redshifts, though this number is also strongly subject to 
cosmic variance effects. 

A simple truncated Gaussian model was fit to the com- 
posite spectrum of our high quality LAE candidates. The 
best-fit velocity dispersion was 136 km s~^, suggesting 
that our galaxies lie at the low end of the observed LAE 
mass distribution. While the model fits reasonably well, 
there were two noticeable discrepancies for which the 
model failed to account. First, we found modest evi- 
dence for excess light at 1214 A and 1215 A, which we 
tentatively attributed to a non-trivial Lya escape frac- 
tion. There was also an observed excess at 1217.5 A that 
may be the result of galactic outflows separated from the 
LAE candidates by 440 km s~^. As this composite rep- 
resents the average properties of our LAE candidates, 
this observed excess implies that outflow processes may 
be prevalent in low mass star-forming galaxies at high 
redshift. 

We find the density of LAEs to be --1.5x10"^ Mpc"^ 
for L(Lya)> 6 x 10'*'^ ergs s~^, a frequency far higher 
than any other search for LAEs at comparable redshifts. 
We find that the excess is instead consistent with extrap- 
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olated cumulative number densities of higher (z > 5.7) 
LAE surveys, initiaUy suggesting minimal evolution in 
the LAE number density between z = 4.8 and z = 5.7. 

We report on the possible discovery of two structures at 
z ~ 4.4 (three members) and z ~ 4.8 (seven members). 
Removing these galaxies from our sample and adopting 
the remaining galaxies as "field" LAEs, we find number 
densities consistent with lower rcdshift (z ^ 4.5) surveys, 
allowing for evolution in the LAE number density. 

We investigate the effect of cosmic variance using 
simulated observations of four samples of narrowband- 
imaging-selectcd LAEs. The results of these simulations 
suggest that wc cannot not rule out cosmic variance as 
the sole cause for the observed excess in the LAE density 
in the C11604 field. The results of the simulations also 
suggest that our field contains a large scale structure of 
LAEs, consistent with the observed redshift clustering of 
the C11604 LAE candidates. 

Best-fit Schechter parameters were determined by sim- 
ulating the effect of observing LAEs with our instrumen- 
tal setup to account for unknown slit attenuation. The 
resultant best-fit Schechter product ($*L*) was found to 
be = 2.2;^:^ X 1039 ergs s"! Mpc'^. Although 

these simulations are sensitive to the assumed complete- 
ness limit of the survey, we find that our results are gen- 
erally consistent with other surveys both at intermediate 
{z ^ 4.5) and high (z ^ 6) redshifts. 

Simulating the weak lensing effect induced by the 
C11604 supercluster using SIS models characterized 
by published cluster velocity dispersions, we find an 
average increase of 6% in the observed number counts 
(or equivalently luminosity) of simulated LAE pop- 
ulations between z = 4.4 and z = 6.5 as compared 
to unlensed populations. The observed change in the 
best-fit product of the luminosity function parameters 
($,L*) due to weak lensing was 7.5%, consistent within 
the errors to the unlensed values, and far too small to 
explain our observed number density excess. We also 
investigated the effects of strong lensing induced by the 
supercluster, finding that no galaxies are likely strongly 
lensed by the cluster potential. 
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Fig. 19. — Spectral ID, cutout of flux calibrated DEIMOS one-dimensional spectrum uncorrected for slit losses in units of /^Jy, and cutout 
of DEIMOS two-dimensional spectrum for each LAE candidate. Postage stamps of the ACS F606W and F814W images (when available) 
or the LFC r' , i' , and z' show the DEIMOS slit (box) and the LAE candidate position (circle) either from the detected position or inferred 
assuming the LAE fell in the middle of the slit (widthwise). LAE candidates 1 through 6. 
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Fig. 20.— LAE candidates 7 through 12. 
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Fig. 21. — LAE candidates 13 through 17. 
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,95; 


SC1NM1.85 


4.82165 


241, 


.073486 


43, 


.310581 


1 


4.92t^|« 


1, 


.18: 


SC1NM2.79 


4.84754 


241, 


.066261 


43, 


.310909 


1 


071+2.17 

'^•'^-1.43 


0, 


.90: 


SC2NM1.34 


5.66885 


241, 


.118866 


43, 


.204311 


3 




1, 


.62; 


SC2NM1.45 


5.62570 


241, 


.196610 


43, 


,212936 


3 


5.56t|?i 


1, 


.91: 


SC2NM1.78 


4.84561 


241, 


.140350 


43, 


.334343 


1 


c 00 + 3,63 
°-^^-2.38 


1, 


.51: 


SC2NM2.61 


4.85259 


241, 


.181778 


43, 


.274860 


3 


^•^^-1.19 


0, 


.72; 


16XR1.26 


4.39288 


241, 


.089600 


43, 


.244692 


3 


3.68li-57 


0, 


.71: 


16XR1.97 


5.02973 


241, 


.094269 


43, 


,292652 


1 


8.97i^j° 


2, 


,37; 


16XR2.19 


4.53375 


240, 


.858124 


43, 


,356789 


3 


2.841J-6J 


0, 


,59; 


16XR3.26 


4.81549 


240, 


.968445 


43, 


.382172 


3 


6.75l^;|? 


1, 


.61: 


16XR3.44 


4.69732 


240, 


.017355 


43. 


.286700 


2 


4.90l|?^ 


1. 


11: 



,+2.18 
'—1,40 


t- „+3.16 


25.24 


24.79 


25, 


.63 


25, 


.34 


24. 


,03 


>28.7 


>143.4 


>19.0 


,+ 1.45 

'-0.93 


4 m+2-30 


25.28 


24.88 


25, 


.72 


25, 


.20 


24. 


.01 


>13.6 


>65.4 


>11.9 


+ 1.12 
'-0.72 


2 q2+^-'^'^ 
^■'^^-i.os 


25.36 


24.83 


25, 


.00 


25, 


.13 


23, 


,39 


>12.1 


>55.1 


>9.5 


,+0.69 
'-0.45 


1 07+1,09 

J^-°'-0.72 


25.40 


24.85 


25, 


.69 


25, 


.32 


23, 


,89 


>6.7 


>27.7 


>4.8 


1+0.53 
-0.35 


1 44+0-84 
^•*^-0.55 


24.98 


24.58 


25, 


.72 


25, 


.31 


23. 


,98 


>4.0 


>18.8 


>3.1 


,+0.95 
'-0.63 


^■"'-'-1.14 


25.41 


24.87 


25, 


.21 


25, 


.31 


23. 


,61 


>4.5 


>21.4 




+ 1.12 
-0.74 


3.45t?;S 


i 


I 


26 


.18 


25 


.22 


23, 


,67 


>6.6 


m 


k 


+0.88 

-0,57 


2.40lJ:« 


26.11 


25.62 


25, 


.59 


25, 


.17 


23, 


,93 


>15.0 


>78.5 


>11.9 


,+0.43 
-0.29 


1 14+0.68 
^■^^-0.46 


I 


I 


25, 


.94 


24, 


.87 


23. 


.65 


>4.6 


m 


>2.1 


+0.44 
-0.30 




25.25 


24.70 


25, 


.60 


25, 


.28 


23. 


,94 


>3.2 


>16.3 


>2.5 


.+ 1.38 

-0,90 




25.44 


24.92 


25, 


.57 


25 


.37 


23, 


,96 


>11.5 


>47.1 


>8.5 


,+0.35 
-0,23 




25.43 


24.94 


25, 


.43 


24, 


.91 


23, 


,86 


>3.3 


>13.7 


>2.5 


+0.94 
-0.61 


2 56+1-49 


25.50 


24.95 


25, 


.60 


25, 


.21 


23, 


,86 


>6.0 


>23.6 


>4.5 


+0.69 
-0.49 


71+1.06 
-I- '-1-0.75 


25.47 


24.88 


25, 


.68 


25, 


.27 


23. 


,94 


>4.3 


>17.0 


>3.2 



IDs Eire generated from a combination of mask names and slit numbers. 
''Computed assuming the LAE candidate falls in the center of the slit (widthwise). 
Confidence in a candidate as a genuine LAE. 3 is most secure, 1 is least secure. 

Lower limit, calculated with a slit throughput of 0.8. Errors include a ~60% systematic uncertainty, which result from uncertainties in absolute flux measurements of DEIMOS spectra. 
Calculated using a line flux corrected for attenuation of Lya photons due to intervening HI regions. 
f All ACS and LFC magnitudes arc 3a- limiting magnitudes. ACS 3o- magnitudes are calculated using a 0.42" circular aperture and LFC 3(T magnitudes are calculated using a 1" circular 
aperture. 

^ Rest-frame EW. Calculated using the 3(t limiting magnitude in a band encompassing the Lya line. The ACS 3<t limiting magnitudes were used when available. 

Rest-frame EW. Calculated using the turnover magnitude (sec §4.1) in a band encompassing the Lya line. 
' Rest-frame EW. Calculated using the 3o- limiting magnitude in a band redward of the Lya line. 
^ Rest-frauic EW. Calculated using the turnover magnitude in a band redward of the Lya line. 
*^ No bands completely redward of the Lya emission. 
' Outside the ACS coverage. 

""Not calculated for LAEs outside the ACS coverage as the turnover magnitude was similar to the 3<t limiting magnitude in LFC images 



00 



TABLE 2 

Properties of LAE candidates detected in the imaging 



ID 


z 


0:2000 


52000 


Class FhyoL'^ L^yct"' 


SFR 


SFR(7v " M[/^^= mFeoew "t-f814w 


r' 


i' z' 


EW ''• 












(10-18 ergs s-i cm-2) (lO*^ erg s'^) 


M© yr-i 


Mo yr-i 




(A) 


(A) 





FG1.20 4.40149 241.066400 43.343734 2 

GHF2.61S.3 1.85025 240.991890 43.413297 2 
16XR1.72 4.71230 241.062630 43.323880 3 



75.42l«;°| 



1 71+0.99 

17.12tl-ll 



2 47+1-43 
^•*'_0.94 



6.24 
26.52 



+3.58 
-2.28 
15.14 
9.65 



14.0 

9.91 



29.75 



+1.43 

1.30 

+0.88 
-0.93 
+0.79 
0.84 



-21.02 

-20.65 
-21.85 



25.53 

26.89 
25.71 



25.16 25.31 
25.70 ..J 
24.45 25.43 



24.41 

..J 
24.11 



... 

..J 
24.17 



>11.0 >12 

>37.3 >66 
>60.6 >152 



" Lower limit, calculated with a slit throughput of 0.8. Errors include a ~60% systematic uncertainty, which result from uncertainties in absolute flux measurements of DEIMOS spectra. 
''Calculated using a line flux corrected for attenuation of Lya photons due to intervening HI regions. 
'^Calculated from the F814W magnitude, roughly a measure of the flux density near rest-frame 1500 A 
''Rest-frame EW. Calculated using the magnitude in a band encompassing the Lya line. 

Rest-frame EW. Calculated using the magnitude in a band redward of the Lya line. 
■^Not detected. 
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